SPECIFICATION 



NON-LINEAR CHARACTERISTIC REPRODUCING APPARARUS AND NON- 
LINEAR CHARACTERISTIC REPRODUCING PROGRAM STORAGE MEDIUM 

5 

BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to a non- linear 
characteristic reproducing apparatus for performing a 

10 simulation in which a non-linear transformation is applied 
to a state quantity and then outputted, and a non- linear 
characteristic reproducing program storage medium storing 
therein a non-linear characteristic reproducing program 
which causes a computer to operate as a non-linear 

15 characteristic reproducing apparatus . 
Description of the Related Art 

r - The great majority of products and components, 

which are an object of design and development, includes a 
non-linear characteristic. This is a subject which cannot 

20 be avoided throughout all the processes of a product 

development from planning via a design to a test of trial 
manufacture. However, the non-linearity, which has an 
effect on all parts of the product development, has two 
aspects one of which is an obstacle to the product 

25 development and another effective, and has very important 
meaning. With respect to the aspect which is an obstacle 
to the product development, it is a cause of the uncertain 
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phenomenon and in many cases it happens that the uncertain 
phenomenon has a bad effect on a function of manufactured 
goods and deteriorates its performance and reliability. 
Thus, such a non-linearity has to be removed or has to 
5 avoid an influence on manufactured goods . Next , with 

respect to the effective aspect of the non-linearity, it 
often happens that the non-linearity is utilized actively 
on purpose as means of creation and implementation of 
functions and characteristics, such as a non-linear spring, 
10 a semiconductor, a link mechanism and clutch mechanism, 

which are not be realized by use of the linearity. In this 
case, application of the non-linearity is indispensable for 
implementation of the function and achievement of the 
performance . 

15 To provide a model, there is a need to reproduce 

faithfully the non-linearity. However, most of the 
conventional techniques for providing a model, for example, 
the finite element method, is originally developed for a 
linear system. Therefore, generally it is difficult and 

20 troublesome to deal with the non-linearity in accordance 
with the conventional techniques for providing a model. 
Consequently, in the conventional product model, it is 
usual that the non-linearity is omitted from beginning, or 
alternatively the non-linearity is dealt with as the 

25 equivalent linearity in which the non-linearity is 

approximated or averaged. In the event that the non- 
linearity greatly effects on the performance, or in the 
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event that the non-linearity is utilized to implement the 
necessary function, hitherto it was obliged to use the 
numerical processing method which is complicated and 
lacking of generality, so alternatively we had to cope with 
5 those cases individually in accordance with intuition and 
experience of a skilled engineer. This makes it difficult 
to provide a model. And also this is a primary factor to 
prevent a computer from being utilized for a product 
development . Indeed an actual product is involved in 
10 various types of complicated non-linearity, and thus it is 
difficult to solve those various types of non-linearity in 
accordance with a common theory. 



SUMMARY OF THE INVENTION 
15 In view of the foregoing, it is an object of the 

present invention to provide a non-linear characteristic 
reproducing apparatus for reproducing a non- linear behavior 
of manufactured goods and components on a simulation, and a 
non- linear characteristic reproducing program storage 
20 medium storing therein a non-linear characteristic 

reproducing program which causes a computer to operate as a 
non-linear characteristic reproducing apparatus. 

To achieve the above-mentioned object, the present 
invention provides first non-linear characteristic 
25 reproducing apparatus wherein upon receipt of a 
predetermined first state quantity, a non- linear 
transformation processing is applied to the entered first 
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state quantity so that a second state quantity is generated 
and outputted, said non-linear characteristic reproducing 
apparatus comprising: 

a state quantity transformation unit for linear- 
5 transforming the first state quantity to the second state 
quantity every sampling time in accordance with a 
transformation parameter set up; and 

a non-linear characteristic reproducing unit for 
determining a transformation parameter for transformation 

10 at a subsequent sampling time in accordance with an 

estimated observation quantity at the subsequent sampling 
time of at least one state quantity of the first state 
quantity and the second state quantity or a state quantity 
derived from said one state quantity to set the determined 

15 transformation parameter to said state quantity 
transformation unit. 

^ . In the first non-linear characteristic reproducing 

apparatus according to the present invention, it is 
acceptable that said non- linear characteristic reproducing 
20 unit receives the estimated observation quantity and one or 
more variables as well, and determines the transformation 
parameter in accordance with the estimated observation 
quantity and one or more variables thus received. 

Alternatively, it is acceptable that said non- 
25 linear characteristic reproducing unit receives the 

estimated observation quantity and one or more variables as 
well, and determines a plurality of second state quantities 
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in accordance with one variable. 

In the first non-linear characteristic reproducing 
apparatus according to the present invention, it is 
preferable that said non-linear characteristic reproducing 
5 unit determines in form of the transformation parameter a 

normalized estimation value in which an estimation value of 
the second state quantity at the subsequent sampling time 
is normalized with the estimated observation quantity. In 
this case, it is preferable that said non- linear 

10 characteristic reproducing unit the estimation value of the 
second state quantity at the subsequent sampling time is 
divided or differentiated with an absolute value of the 
estimated observation quantity, so that the estimated 
observation quantity is determined in form of the 

15 transformation parameter. 

In the first non-linear characteristic reproducing 
apparatus according to the present invention, it is 
preferable that said non-linear characteristic reproducing 
apparatus is an apparatus for reproducing characteristics 

20 of a non -linear spring, 

said state quantity transformation unit performs a 
transformation between a velocity difference on both ends 
of the non -linear spring and a variation of load of the 
non -linear spring, and 

25 said non-linear characteristic reproducing unit 

determines a transformation parameter for performing 
transformation between the velocity difference on both ends 
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of the non-linear spring as an object of a characteristic 
reproduction and the variation of load of the non- linear 
spring at a subsequent sampling time in accordance with an 
estimated observation quantity at the subsequent sampling 
5 time of the velocity difference on both ends of the non- 
linear spring to set the determined transformation 
parameter to said state quantity transformation unit. 

Alternatively, in the first non- linear 
characteristic reproducing apparatus according to the 
10 present invention, it is ' preferable that said non-linear 
characteristic reproducing apparatus is an apparatus for 
reproducing characteristics of an air spring, 

said state quantity transformation unit performs a 
transformation between a deformation velocity of the air 
15 spring and a variation of an internal pressure of the air 
spring, and 

„ . said non-linear characteristic reproducing unit 

determines a transformation parameter for performing 
transformation between the deformation velocity of the air 

20 spring as an object of a characteristic reproduction and 

the variation of the internal pressure of the air spring at 
a subsequent sampling time in accordance with an estimated 
observation quantity at the subsequent sampling time of the 
deformation velocity of the air spring to set the 

25 determined transformation parameter to said state quantity 
transformation unit. 

Still alternatively, in the first non-linear 
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characteristic reproducing apparatus according to the 
present invention, it is preferable that said non-linear 
characteristic reproducing apparatus is an apparatus for 
reproducing characteristics of a link mechanism, 
5 -said state quantity transformation unit transforms 

a value of velocity or angular velocity of a supporting 
point portion of the link mechanism, and 

said non- linear characteristic reproducing unit 
determines a transformation parameter for transforming the 

10 value of velocity or angular velocity of the supporting 
point portion of the link mechanism as an object of a 
characteristic reproduction at a subsequent sampling time 
in accordance with an estimated observation quantity at the 
subsequent sampling time of velocity or angular velocity 

15 applied to the supporting point portion of the link 

mechanism to set the determined transformation parameter to 
said state quantity transformation unit. 

In this case, it is acceptable that said state 
quantity transformation unit transforms a value of velocity 

20 or angular velocity of the supporting point portion of the 
link mechanism, and also a value of force or torque applied 
to the supporting point portion, using the same parameter 
as that used for transformation of the value of velocity or 
angular velocity of the supporting point portion of the 

25 link mechanism, said parameter being set by said non-linear 
characteristic reproducing unit. 

Further in the first non-linear characteristic 
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reproducing apparatus according to the present invention, 
it is preferable that said non-linear characteristic 
reproducing apparatus is an apparatus for reproducing 
characteristics of an object moving while involving a 
5 friction, 

said state quantity transformation unit performs a 
transformation between a force applied to the object and a 
moving velocity of the object, and 

said non- linear characteristic reproducing unit 

10 determines a transformation parameter for performing 

transformation between the force applied to the object as 
an object of a characteristic reproduction and the moving 
velocity of the object at a subsequent sampling time in 
accordance with an estimated observation quantity at the 

15 subsequent sampling time of the moving velocity of the 

object to set the determined transformation parameter to 
said state quantity transformation unit . 

In this case, it is acceptable that said non- 
linear characteristic reproducing unit determines a 

20 frictional force applied to the object in form of the 
transformation parameter to set the determined 
transformation parameter to said state quantity 
transformation unit, and 

said state quantity transformation unit performs a 

25 transformation between the force applied to the object and 
the moving velocity of the object, in a case where a force, 
in which the frictional force set by said state quantity 



transformation unit is subtracted from the force applied to 
the object, is applied to an object in which an effect of 
the frictional force is neglected. 

In this case, it is acceptable that said non- 
5 linear characteristic transformation unit comprises: a 

kinetic friction generating unit for determining a kinetic 
frictional force applied to the object as an object of a 
characteristic reproduction at a subsequent sampling time 
in accordance with an estimated observation quantity at the 

10 subsequent sampling time' of the moving velocity of the 

object; a static friction generating unit for determining a 
static frictional force applied to the object as an object 
of a characteristic reproduction at a subsequent sampling 
time in accordance with an estimated observation quantity 

15 at the subsequent sampling time of the force applied to the 
object; and a frictional force selection unit for selecting 
one frictional force between the dynamic frictional force 
generated in said kinetic friction generating unit and the 
static frictional force generated in said static friction 

20 generating unit in accordance with the estimated 

observation quantity at the subsequent sampling time of the 
moving velocity of the object to set the selected 
frictional force to said state quantity transformation unit. 
Further, in this case, it is acceptable that said 

25 kinetic friction generating unit determines a kinetic 

frictional force different in value in accordance with an 
estimated observation quantity at the subsequent sampling 
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time of the moving velocity of the object. 

In the first non- linear characteristic reproducing 
apparatus according to the present invention, it is 
preferable that said non-linear characteristic reproducing 
5 apparatus is an apparatus for reproducing characteristics 
of a variable inertia moment mechanism in which an inertia 
moment is varied in accordance with an angular velocity, 

said state quantity transformation unit performs a 
transformation between a torque applied to the variable 

10 inertia moment mechanism' and an angular velocity of the 
variable inertia moment mechanism, and 

said non- linear characteristic reproducing unit 
determines a transformation parameter for performing 
transformation between the torque applied to the variable 

15 inertia moment mechanism as an object of a characteristic 
reproduction and an angular acceleration velocity of the 
variable inertia moment mechanism at a subsequent sampling 
time in accordance with an estimated observation quantity 
at the subsequent sampling time of the angular velocity of 

20 the variable inertia moment mechanism to set the determined 
transformation parameter to said state quantity 
transformation unit. 

In this case, it is acceptable that said variable 
inertia moment mechanism as an object of a characteristic 

25 reproduction has a translational member translating in a 
radius direction in accordance with a centrifugal force, 
and 
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said non- linear characteristic reproducing unit 
comprises: a rotational translation transformation unit for 
determining a centrifugal force applied to the 
translational member in accordance with an estimated 
5 observation quantity at a subsequent sampling time of the 
angular velocity of the variable inertia moment mechanism; 
and a translational motion reproducing unit for reproducing 
a translational motion of the translational member by the 
centrifugal force determined by said rotational translation 
10 transformation unit, and' said non-linear characteristic 
reproducing unit determines a transformation parameter 
according to the translational motion of the translational 
member reproduced by said translational motion reproducing 
unit to set the determined transformation parameter to said 
15 state quantity transformation unit. 

To achieve the above-mentioned object, the present 
invention provides a second non- linear characteristic 
reproducing apparatus comprising a linear model unit for 
reproducing characteristics of a linear system, and a non- 
20 linear model unit for determining, upon receipt of an 

estimated observation quantity at a subsequent sampling 
time of at least one state quantity of a first state 
quantity and a second state quantity, which are in a 
relation of mutually non-linear transformation, or a state 
25 quantity derived from said one state quantity, from said 
linear model unit, a transformation parameter used for a 
linear transformation at the subsequent sampling time 
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between the first state quantity and the second state 
quantity, 

wherein an operation of determining the 
transformation parameter at the subsequent sampling time in 
5 said non-linear model unit and a linear operation including 
a linear transformation between the first state quantity 
and the second state quantity using the transformation 
parameter at the subsequent sampling time determined in 
said non-linear model unit, in said linear model unit are 

10 alternatively repeated. 1 

In the second non- linear characteristic 
reproducing apparatus according to the present invention, 
it is acceptable that said non-linear model unit has a 
plurality of non-linear transformation units for 

15 determining, upon receipt of an estimated observation 
quantity at a subsequent sampling time of at least one 
state quantity of a first state quantity and a second state 
quantity, which are in a relation of mutually non -linear 
transformation and at least one of which is different in 

20 type, or a state quantity derived from said one state 
quantity, a transformation parameter used for a linear 
transformation at the subsequent sampling time between the 
first state quantity and the second state quantity. 
In the second non-linear characteristic 

25 reproducing apparatus according to the present invention, 
it is preferable that said non -linear characteristic 
reproducing apparatus is an apparatus for reproducing 
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characteristics of Geneva mechanism, 

said linear model unit performs linear 
transformations including a linear transformation of 
angular velocity- to-angular velocity between master section 
5 and slave section of the Geneva mechanism, and a linear 

transformation of torque-to-torque between master section 
and slave section of the Geneva mechanism, and 

said non-linear model unit has a non-linear 
transformation unit for determining a transformation 

10 parameter including information as to whether master 
section and slave section of the Geneva mechanism are 
connected to one another at the subsequent sampling time, 
said information being used for both the linear 
transformation of angular velocity- to-angular velocity 

15 between master section and slave section of the Geneva 

mechanism and the linear transformation of torque-to-torque 
between master section and slave section of the Geneva 
mechanism, at a subsequent sampling time of an angle of the 
master section of the Geneva mechanism, in accordance with 

20 an estimated observation quantity at the subsequent 
sampling time. 

In the second non-linear characteristic 
reproducing apparatus according to the present invention, 
it is preferable that said non- linear characteristic 

25 reproducing apparatus is an apparatus for reproducing 

characteristics of a liquid residue warning lamp in which a 
lamp and a thermistor are connected in series. 
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said linear model unit performs linear 
transformations including a transformation between a 
voltage applied to the lamp and a current conducting 
through the lamp and a transformation between a voltage 
5 applied to the thermistor and a current conducting through 
the thermistor, and 

said non-linear model unit has a first non-linear 
transformation unit for determining a resistance of the 
thermistor at a subsequent sampling time of a consumed 
10 power of the thermistor in accordance with an estimated 
observation quantity at the subsequent sampling time. 

To achieve the above-mentioned object, the present 
invention provides a third non-linear characteristic 
reproducing apparatus comprising: 
15 a logical decision unit for receiving one or more 

variables , and determining a logical value at a subsequent 
sampling time, selected among from a plurality of discrete 
values in accordance with the received one or more 
variables ; and 

20 a state quantity selecting unit for receiving a 

predetermined input state quantity and outputting at the 
subsequent sampling time an output state quantity in which 
a relation between the input state quantity and the output 
state quantity is changed over to a relation according to 

25 the logical value at the subsequent sampling time 
determined by said logical decision unit . 

Here, the above-mentioned "plurality of discrete 
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values" may be, for example, binary ('0' and ' 1') or 
ternary ('0' and " 1 ' and '-1'). 

Further, the above-mentioned "the first state 
quantity is changed over" implies that the first state 
5 quantity is changed over to various states, for example, a 
selection between passage and block of the first state 
quantity in the event that the logical value is expressed 
by the binary, and for example, a selection among forward 
flow, reverse flow and flow stop of the first state 

10 quantity in the event that the logical value is expressed 
by the ternary. 

In a third non-linear characteristic reproducing 
apparatus according to the present invention, it is 
acceptable that said state quantity selecting unit outputs 

15 at the subsequent sampling time an output state quantity in 
which a relation between the input state quantity and the 
output state quantity is changed over to a connection 
relation according to the logical value at the subsequent 
sampling time determined by said logical decision unit. 

20 Alternatively it is acceptable that said state 

quantity selecting unit integrates the input state quantity 
and outputs as an output state quantity, and in a case 
where the logical value at the subsequent sampling time 
determined by said logical decision unit is a predetermined 

25 logical value, said state quantity selecting unit outputs 
an output state quantity which is changed over to an 
initial value at the subsequent sampling time. 
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Further, it is acceptable that said logical 
decision unit determines a logical value at a subsequent 
sampling time, selected among from a plurality of discrete 
values in accordance with an estimated observation quantity 
5 at the subsequent sampling time of a plurality of input 
state quantities, and 

said state quantity selecting unit receives the 
plurality of input state quantities and outputs at the 
subsequent sampling time in form of an output state 
10 quantity an input state quantity selected in accordance 
with the logical value at the subsequent sampling time 
determined by said logical decision unit. 

In the third non- linear characteristic reproducing 
apparatus according to the present invention, it is 
15 preferable that said non-linear characteristic reproducing 
apparatus is an apparatus for reproducing characteristics 
of- a mechanism including two members having a relative 
movement possible state and a relative movement impossible 
state, 

20 said logical decision unit determines a logical 

value representing whether said two members are in the 
relative movement possible state or the relative movement 
impossible state in accordance with an estimated 
observation quantity at a subsequent sampling time of a 

25 relative position of said two members or a state quantity 

derived from the relative position of said two members, and 
said state quantity selecting unit changes over a 



16 



relation between a relative moving velocity and a 
distribution of force of said two members . 

In the third non-linear characteristic reproducing 
apparatus according to the present invention, it is 
5 preferable that said non-linear characteristic reproducing 
apparatus is an apparatus for reproducing characteristics 
of a clutch mechanism including two members having a 
relatively slid sliding state and a mutually connected 
connecting state, 

10 said logical decision unit determines a logical 

value representing whether said two members are in the 
sliding state or the connecting state in accordance with an 
estimated observation quantity at a subsequent sampling 
time of a relative angular velocity of said two members, 

15 and 

said state quantity selecting unit changes over a 
relation between a relative angular velocity and a 
distribution of torque of said two members . 

Further, in the third non- linear characteristic 

20 reproducing apparatus according to the present invention, 
it is preferable that said non -linear characteristic 
reproducing apparatus is an apparatus for reproducing 
characteristics of a brake mechanism for applying a braking 
energy to a driving shaft, 

25 said logical decision unit determines a logical 

value representing whether the driving shaft is in a 
rotating state or a stationary state in accordance with an 
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estimated observation quantity at a subsequent sampling 
time of a torque which a brake receives from the driving 
shaft, and 

said state quantity selecting unit changes over a 
5 braking torque to be applied to the driving shaft . 

In the third non- linear characteristic reproducing 
apparatus according to the present invention, it is 
preferable that said non-linear characteristic reproducing 
apparatus is an apparatus for reproducing characteristics 
10 of an automatic-reset mechanism in which a spring is 

effected on a movable member, and when application of an 
external force to the movable member is removed, the 
movable member is automatically reset to an initial state 
by an effect of the spring, said automatic-reset mechanism 
15. having a stopper to limit a movable range of the movable 
member , 

^ _ said logical decision unit determines a logical 

value representing whether the movable member interferes 
with the stopper in accordance with an estimated 
20 observation quantity at a subsequent sampling time of a 
moving position of the movable member, and 

said state quantity selecting unit changes over a 
relation between velocity or angular velocity of the 
movable member and force or torque to be applied to the 
25 movable member. 

To achieve the above-mentioned object, the present 
invention provides fourth non-linear characteristic 
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reproducing apparatus wherein upon receipt of a 
predetermined first state quantity, a second state quantity 
changed over to a non-linearity according as the received 
first state quantity is in a predetermined state is 
5 generated, said non- linear characteristic reproducing 
apparatus comprising: 

a state variation estimation unit for predicting a 
state quantity variation width during a period from a 
present sampling time to a subsequent sampling time of the 

10 first state quantity; 

a state deviation detection unit for determining a 
deviation between a value at the present sampling time of 
the first state quantity and a decision value for deciding 
whether the first state quantity is in a predetermined 

15 state; 

a stable state decision unit for determining a 
locjical value to be selected from among a plurality of 
discrete values, predicting a non- linear variation at a 
subsequent sampling time in accordance with a comparison of 
20 the state quantity variation width predicted by said state 
variation estimation unit with the deviation determined in 
said state deviation detection unit; and 

a state quantity selecting unit for outputting the 
second state quantity changed over in accordance with the 
25 logical value determined in said stable state decision unit. 

In the fourth non- linear characteristic 
reproducing apparatus according to the present invention. 
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it is acceptable that said state deviation detection unit 
determines a deviation between a value at the present 
sampling time of the first state quantity and a decision 
value for deciding whether the first state quantity is in a 
5 predetermined state, said decision value being varied in 

accordance with positive or negative of the state quantity 
variation width predicted in said state variation 
estimation unit . 

To achieve the above-mentioned object, the present 

10 invention provides a fifth non-linear characteristic 
reproducing apparatus comprising: 

a linear model unit for reproducing 
characteristics of a linear system, including a state 
quantity transformation unit for linear- transforming a 

15 first state quantity to a second state quantity in 

accordance with a transformation parameter set up; and 
r . a non- linear model unit for generating, upon 

receipt of an estimated observation quantity at a 
subsequent sampling time of a predetermined first 

20 observation state quantity from said linear model unit, the 
transformation parameter in accordance with the received 
estimated observation quantity and setting the generated 
transformation parameter on said state quantity 
transformation unit, 

25 wherein said non-linear model unit comprises: 

a slow change reproducing unit for receiving from 
said linear model an observation quantity or an estimated 
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observation value of a predetermined second observation 
state quantity identical to or different from the first 
observation state quantity, to generate a slow change state 
quantity reflecting characteristic of a first non-linear 
5 system offering a relatively slow behavior change in 

accordance with the observation quantity or the estimated 
observation value of the second observation state quantity; 
and 

a characteristic generating unit for generating a 
10 transformation parameter- reflecting characteristic of a 
second non-linear system offering a relatively rapid 
behavior change in accordance with the estimated 
observation value at a subsequent sampling time of the 
first observation state quantity derived from said linear 
15 model unit and the slow change state quantity generated in 
said slow change reproducing unit and setting up the 
generated transformation parameter on said state quantity 
transformation unit . 

In this case, it is preferable that said slow 
20 change reproducing unit comprises: 

a stationary value setting up unit for determining 
a stationary value, after passage of an infinite time, of 
the slow change state quantity assuming that a state of 
said linear model unit is maintained, in accordance with 
25 the observation quantity or the estimated observation value 
of the second observation state quantity; and 

a normalization response unit for generating the 
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slow change state quantity reflecting characteristic of the 
first non-linear system, which is to be transferred to said 
characteristic generating unit, in accordance with a 
stationary value of the slow change state quantity 
5 determined in said stationary value setting up unit, and a 
known normalized time change characteristic of the slow 
change' state quantity in the first non- linear system. 

in this case, it is preferable that said non- 
linear characteristic reproducing apparatus is an apparatus 

10 for reproducing characteristics of a system having an 

element which is varied in a resistance value in accordance 
with a temperature variation, 

said state quantity transformation unit sets up a 
transformation parameter representative of the resistance 

15 value and performs a transformation between a voltage 

applied to the element having the resistance value and a 
current conducting through the element, 

said stationary value setting up unit determines 
a stationary value, after passage of an infinite time, of a 

20 temperature of the element assuming that an energy to be 
consumed in the element is maintained, in accordance with 
the estimated observation quantity at a subsequent sampling 
time of the energy to be consumed in the element, 

said normalization response unit determines a 

25 temperature of the element in accordance with a stationary 
value of the temperature of the element determined in said 
stationary value setting up unit, and a known normalized 
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step response curve representative of a time variation of 
the temperature of the element to a variation of an energy 
to be consumed in the element, and 

said characteristic generating unit generates a 
5 transformation parameter representative of the resistance 
value of the element in accordance with the temperature 
determined in said normalization response unit and sets up 
the generated transformation parameter on said state 
quantity transformation unit. 

10 In this case, it is acceptable that said 

normalization response unit determines a temperature of the 
element in accordance with an observation quantity of an 
ambient temperature . 

Incidentally, for example, velocity and angular 

15 velocity or force and torque are determined according as 

the mechanism and the like as an object of reproduction is 
a translational system or a rotational system. 
Consequently, in any of the present invention as mentioned 
above, even if it is not clearly stated, force is a concept 

20 including torque, torque is a concept including force, 

velocity is a concept including angular velocity, angular 
velocity is a concept including velocity, position is a 
concept including rotational angle. Other physical 
quantity is also interpreted in a similar fashion. 

25 To achieve the above-mentioned object, the present 

invention provides a first non- linear characteristic 
reproducing program storage medium storing a non-linear 
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characteristic reproducing program which causes a computer 
to operate as a non- linear characteristic reproducing 
apparatus wherein upon receipt of a predetermined first 
state quantity, a non-linear transformation processing is 
5 applied to the entered first state quantity so that a 
second state quantity is generated and outputted, 

wherein said non-linear characteristic reproducing 
program storage medium stores a non-linear characteristic 
reproducing program comprising: 

10 a state quantity transformation unit for linear- 

transforming the first state quantity to the second state 
quantity every sampling time in accordance with a 
transformation parameter set up; and 

a non- linear characteristic reproducing unit for 

15 determining a transformation parameter for transformation 
at a subsequent sampling time in accordance with an 
estimated observation quantity at the subsequent sampling 
time of at least one state quantity of the first state 
quantity and the second state quantity or a state quantity 

20 derived from said one state quantity to set the determined 
transformation parameter to said state quantity 
transformation unit . 

To achieve the above-mentioned object, the present 
invention provides a second non-linear characteristic 

25 reproducing program storage medium storing a non-linear 

characteristic reproducing program which causes a computer 
to operate as a non- linear characteristic reproducing 
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apparatus for reproducing characteristics of a system 
including a non- linear system, 

wherein said non- linear characteristic reproducing 
program storage medium stores a non-linear characteristic 
5 reproducing program comprising: 

a linear model unit for reproducing 
characteristics of a linear system, and a non-linear model 
unit for determining, upon receipt of an estimated 
observation quantity at a subsequent sampling time of at 

10 least one state quantity of a first state quantity and a 

second state quantity, which are in a relation of mutually 
non- linear transformation, or a state quantity derived from 
said one state quantity, from said linear model unit, a 
transformation parameter used for a linear transformation 

15 at the subsequent sampling time between the first state 
quantity and the second state quantity, 
r . wherein an operation of determining the 

transformation parameter at the subsequent sampling time in 
said non- linear model unit and a linear operation including 

20 a linear transformation between the first state quantity 
and the second state quantity using the transformation 
parameter at the subsequent sampling time determined in 
said non-linear model unit, in said linear model unit are 
alternatively repeated. 

25 To achieve the above-mentioned object, the present 

invention provides a third non- linear characteristic 
reproducing program storage medium storing a non- linear 
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characteristic reproducing program which causes a computer 
to operate as a non-linear characteristic reproducing 
apparatus for reproducing characteristics of a system 
including a non-linear system, 
5 wherein said non- linear characteristic reproducing 

program storage medium stores a non -linear characteristic 
reproducing program comprising: 

a logical decision unit for receiving one or more 
variables , and determining a logical value at a subsequent 

10 sampling time, selected among from a plurality of discrete 
values in accordance with the received one or more 
variables; and 

a state quantity selecting unit for receiving a 
predetermined input state quantity and outputting at the 

15 subsequent sampling time an output state quantity in which 
a relation between the input state quantity and the output 
state quantity is changed over to a relation according to 
the logical value at the subsequent sampling time 
determined by said logical decision unit. 

20 To achieve the above-mentioned object, the present 

invention provides a fourth non- linear characteristic 
reproducing program storage medium storing a non-linear 
characteristic reproducing program which causes a computer 
to operate as a non- linear characteristic reproducing 

25 apparatus wherein upon receipt of a predetermined first 

state quantity, a second state quantity changed over to a 
non-linearity according as the received first state 
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quantity is in a predetermined state is generated, 

wherein said non- linear characteristic reproducing 

program storage medium stores a non -linear characteristic 

reproducing program comprising: 
5 a state variation estimation unit for predicting a 

state quantity variation width during a period from a 

present sampling time to a subsequent sampling time of the 

first state quantity; 

a state deviation detection unit for determining a 
10 deviation between a value at the present sampling time of 

the first state quantity and a decision value for deciding 

whether the first state quantity is in a predetermined 

state; 

a stable state decision unit for determining a 
15 logical value to be selected from among a plurality of 
discrete values, predicting a non-linear variation at a 
subsequent sampling time in accordance with a comparison of 
the state quantity variation width predicted by said state 
variation estimation unit with the deviation determined in 
20 said state deviation detection unit; and 

a state quantity selecting unit for outputting the 
second state quantity changed over in accordance with the 
logical value determined in said stable state decision unit. 

To achieve the above-mentioned object, the present 
25 invention provides a fifth non-linear characteristic 

reproducing program storage medium storing a non-linear 
characteristic reproducing program which causes a computer 
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to operate as a non- linear characteristic reproducing 
apparatus for reproducing characteristics of a system 
including a non -linear system, 

wherein said non-linear characteristic reproducing 
5 program storage medium stores a non- linear characteristic 
reproducing program comprising: 

a linear model unit for reproducing 
characteristics of a linear system, including a state 
quantity transformation unit for linear- transforming a 

10 first state quantity to a second state quantity in 

accordance with a transformation parameter set up; and 

a non-linear model unit for generating, upon 
receipt of an estimated observation quantity at a 
subsequent sampling time of a predetermined first 

15 observation state quantity from said linear model unit, the 
transformation parameter in accordance with the received 
es-timated observation quantity and setting the generated 
transformation parameter on said state quantity 
transformation unit, 

20 wherein said non-linear model unit comprises: 

a slow change reproducing unit for receiving from 
said linear model an observation quantity or an estimated 
observation value of a predetermined second observation 
state quantity identical to or different from the first 

25 observation state quantity, to generate a slow change state 
quantity reflecting characteristic of a first non-linear 
system offering a relatively slow behavior change in 
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accordance with the observation quantity or the estimated 
observation value of the second observation state quantity; 
and 

a characteristic generating unit for generating a 
5 transformation parameter reflecting characteristic of a 
second non- linear system offering a relatively rapid 
behavior change in accordance with the estimated 
observation value at a subsequent sampling time of the 
first observation state quantity derived from said linear 

10 model unit and the slow change state quantity generated in 
said slow change reproducing unit and setting up the 
generated transformation parameter on said state quantity 
transformation unit. 

Incidentally, while the first to fifth non- linear 

15 characteristic reproducing program storage media of the 

present invention are described above on the basic style, 
the non-linear characteristic reproducing programs stored 
in the non -linear characteristic reproducing program 
storage medium according to the present invention includes 

20 all aspects implementing all types of first to fifth non- 
linear characteristic reproducing program storage media of 
the present invention. 

As mentioned above, according to the present 
invention, it is possible to reproduce a non-linear 

25 behavior of products and parts through modeling. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a view showing functional and 
mechanistic models of a running resistance. 

Fig. 2 is a view showing functional and 
5 mechanistic models of a kinetic frictional force and a 
running resistance. 

Fig. 3 is a view showing a functional model of a 
running resistance. 

Fig. 4 is a view showing a functional model in 
10 which a hierarchical structure is expressed by a relation 
of inclusion. 

Fig. 5 is a perspective view of a computer in 
which a non- linear characteristic reproducing apparatus 
according to one embodiment of the present invention is 
15 implemented. 

Fig. 6 is a schematic diagram showing a hardware 
structure of the computer shown in Fig. 5. 

Fig. 7 is a typical illustration showing a 
structure of a non- linear characteristic reproducing 
20 program stored in a storage medium. 

Fig. 8 is a typical illustration showing a 
structure of a non- linear characteristic reproducing 
program stored in a storage medium. 

Fig. 9 is a typical illustration showing a 
25 structure of a non-linear characteristic reproducing 
program stored in a storage medium. 

Fig. 10 is a typical illustration showing a 
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structure of a non- linear characteristic reproducing 
program stored in a storage medium. 

Fig. 11 is a typical illustration showing a 
structure of a non -linear characteristic reproducing 
5 program stored in a storage medium. 

Fig. 12 is a principle explanatory view of a first 
non- linear characteristic reproducing apparatus according 
to an embodiment of the present invention. 

Fig. 13 is a view showing a non-linear parameter 
10 of a flow storage element . 

Fig. 14 is a view showing a structure of a conical 
spring and a load characteristic. 

Fig. 15 is a view showing a functional model of a 
conical spring. 

15 Fig. 16 is a view showing a structural model of a 

pneumatic cylinder. 

Fig. 17 is a view showing a functional model of a 
pneumatic cylinder. 

Fig. 18 is a view showing a simulation result of a 
20 pneumatic cylinder. 

Fig. 19 is a view showing a P-V * v-V diagram of a 
pneumatic cylinder. 

Fig. 20 is a view showing a structural model of an 
oscillating mechanism. 
25 Fig. 21 is a view showing a functional model of an 

oscillating mechanism. 

Fig. 22 is a view showing a functional model of an 
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oscillating mechanism. 

Fig. 23 is a view showing a simulation result of 
an oscillating mechanism. 

Fig. 24 is a view showing an angular 
5 characteristic of a non-linear transfer factor <& b . 

Fig. 25 is a view showing a structural model of 
friction . 

Fig. 26 is a view showing a functional model of 
friction . 

10 Fig. 27 is a view showing a simulation result of 

friction . 

Fig. 28 is a view showing a structural model of 
the variable moment of inertia. 

Fig. 29 is a view showing functional and 
15 mechanistic models of the variable moment of inertia. 

Fig. 30 is a view showing a simulation result of 
the variable moment of inertia. 

Fig. 31 is a principle explanatory view of a 
second non-linear characteristic reproducing apparatus 
20 according to an embodiment of the present invention. 

Fig. 32 is a view showing a structural model of 
the Geneva mechanism. 

Fig. 33 is a view showing functional and 
mechanistic models of the Geneva mechanism. 
25 Fig. 34 is a view showing a simulation result (1) 

of the Geneva mechanism. 

Fig. 35 is a view showing a simulation result (2) 
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of the Geneva mechanism. 

Fig. 36 is an illustration showing a constitution 
of a residue warning lamp. 

Fig. 37 is a circuit diagram of an electric 
5 circuit of a residue warning lamp. 

Fig. 38 is a view showing characteristics of a 
contactless switch to a resistance variation. 

Fig. 39 is a view showing a functional model of a 
residue warning lamp. 
10 Fig. 40 is a view showing a functional model of a 

residue warning lamp. 

Fig. 41 is a view showing voltage -current 
characteristic model of a warning lamp. 

Fig. 42 is a perspective view of a thermistor .- 
15 Fig. 43 is a view showing a mechanistic model of a 

thermistor . 

«. _ Fig. 44 is a view showing step response 

characteristics of a residue warning lamp. 

Fig. 45 is a view showing response characteristics 
20 of a liquid level oscillation. 

Fig. 46 is a principle explanatory view of a third 
non- linear characteristic reproducing apparatus according 
to an embodiment of the present invention. 

Fig. 47 is a view showing a structural model of a 
25 back- lash. 

Fig. 48 is a view showing functional and 
mechanistic models of a back- lash. 
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Fig. 49 is a view showing a structural model of an 
impact absorption damper. 

Fig. 50 is a view showing functional and 
mechanistic models of an impact absorption damper. 
5 Fig. 51 is a view showing a simulation result of 

an impact absorption damper. 

Fig. 52 is a view showing a structural model of a 

clutch . 

Fig. 53 is a view showing torque transfer 
10 characteristics of a clutch. . 

Fig. 54 is a view showing functional and 
mechanistic models of a clutch. 

Fig. 5 5 is a view showing a clutch model for a 
simulation. 

15 Fig. 56 is a view showing a simulation result of a 

clutch. 

„ _ Fig. 57 is a view showing a structural model of a 

brake . 

Fig. 58 is a view showing a functional model of a 

20 brake. 

Fig. 59 is a view showing a simulation result of a 

brake . 

Fig. 60 is a view showing a construction of an 
automatic -reset mechanism. 
25 Fig. 61 is a view showing a functional model of an 

automatic-reset mechanism. 

Fig. 62 is a view showing non-linear 



34 



characteristics of a damping resistance factor. 

Fig. 6 3 is a view showing a construction of a 
rattle device. 

Fig. 64 is a view showing a functional model of a 
5 rattle device. 

Fig. 6 5 is a block diagram of a fourth non- linear 
characteristic reproducing apparatus according to an 
embodiment of the present invention. 

Fig. 66 is a view useful for understanding a 
10 determination state of the non-linear characteristic 
reproducing apparatus shown in Fig. 65. 

Fig. 67 is a view useful for understanding a 
determination state of the non-linear characteristic 
reproducing apparatus shown in Fig. 65, wherein a 
15 hysteresis is added. 

Fig. 68 is a principle explanatory view of a fifth 
nan- linear characteristic reproducing apparatus according 
to an embodiment of the present invention. 

Fig. 69 is a view showing a functional model of a 
20 coil temperature rise. 

Fig. 70 is a view showing a simulation result of a 
coil temperature rise. 

Fig. 71 is a perspective view of a positive 
characteristic thermistor. 
25 Fig. 72 is a view showing a functional model of an 

overload protection. 

Fig. 73 is a view showing resistance ratio 
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characteristic of a positive characteristic thermistor. 

Fig. 74 is a view showing a functional model of a 
positive characteristic thermistor. 

Fig. 75 is a view showing a functional model of a 
5 motor lock. 

Fig. 76 is a view showing a simulation result of a 
motor lock. 

Fig. 77 is a view showing an execution sequence of 
functional and mechanistic models . 
10 Fig. 78 is a view showing main symbols for 

functional model. 

Fig. 79 is a view showing symbols for non- 
linearity. 

Fig. 80 is a view showing initialization of ah 
15 integration amount (symbol). 

Fig. 81 is a view showing an ON-OFF switch. 
- - Fig. 82 is a view showing a model of a logical 

product . 

Fig. 83 is a view showing a model of a logical 
20 summation. 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
Hereinafter, there will be described non-linearity 
25 and the basic concept of the present invention, and then 
embodiment of the present invention will be described. 

A non-linearity is an extremely general property 
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in such a degree that all systems, machines and structures, 
which are utilized all over the world, include the non- 
linearity. 

For example, as the non-linearity relating to 
5 characteristics and state of quantity models, there may be 
concerned an existing place, a generating cause and a 
dependence factor. Phenomenon and behavior of the non- 
linearity are associated with all sorts and conditions , and 
are difficult to be grasped and to be dealt with in its 
10 entirety. 

® With respect to the existing place, the 
non-linearity exists in many places, for example, inside 
material, coupling section and junction section of parts, 
an engagement, a bearing, a gear, an electromagnet, a 
15 semiconductor, etc. 

(D With respect to the generating cause, it is 
classified into two parts of substances and mediums, and of 
mechanisms and structures. The former includes non-linear 
characteristics such as rigidity, damping, electric 
20 resistance and magnetism, and the later includes a 

structural non-linear characteristic such as a universal 
joint, clearance, backlash and friction. 

(3) With respect to the dependence factor, this 
factor may be classified to one, for example, a running 
25 resistance and an air spring, which depends on a state 

quantity such as a displacement, a velocity, a current and 
a voltage, and one which depends on an environment such as 
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a temperature and a humidity, one which depends on a 
geometric structure such as an oscillating mechanis , and 
one which depends on a structural separation and coupling 
such as a backlash and a switch. 
5 With respect to these various types of non- 

linearity, hitherto there has been introduced and 
explicated many non-linear theories and models. These 
present invention is to provide a system of applying the 
explicated non- linear theories and models to products and 

10 components. 

A model of a product or a component , which include 
a non- linear system, comprises of course linear and non- 
linear models. Of the linear and non -linear models, the 
linear mode is implemented in accordance with a linear 

15 theory. In this case, it is a condition that the linear 
mode can be expressed by a mathematical model capable of 
being handled by a lineax algebra. On the other hand, with 
respect to the non- linear function, upon taking these 
linear properties into consideration, there is a need to 

20 implement a model from the following two aspects . 

One of the two aspects is to implement a model 
concerned with the non -linear behavior, in which a 
characteristic, a factor and an attached load (parameter in 
a lump) are affected by state quantity and change from 

25 linear constants in processes wherein features and details 
of mechanisms and physical phenomena of products and 
components are clarified in form of their functions . This 
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non-linearity is affected by also an influence of an 
external environment such as atmospheric pressure and a 
temperature by which the system is suffered. Here such a 
non-linearity is referred to as a non-linear parameter. 
5 Another of the two aspects is a non-linearity 

caused by structures and mechanisms of products and 
components, for example, a structural non-linearity in 
which a compression spring becomes a rigid-body with a 
bottom, and a structural non-linearity in which when a 
10 clutch for a mobile car is joining, an engine is integrated 
with a car body, and when the clutch is opening, both the 
engine and the car body moves separately each other. This 
non-linearity is a model of non-linearity in which the 
parameter and the model itself are changed in such a manner 
15 that the state quantity inside the model is directly 

connected or disconnected. Here such a non-linearity is 
referred to as a structural non-linearity. 

Thus, from the above, the non- linear models can be 
classified to the non-linear parameter in which a parameter 
20 value is varied in accordance with the state quantity, and 
the structural non-linearity in which a model structure is 
varied in accordance with connection and disconnection of 
the state quantity. 

(Technique of implementing a model of non- 
25 linearity) 

( 1 ) A mechanism model 

In order to implement a model of non-linearity, it 
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is possible to implement non-linearity of the model in such 
a manner that a mechanistic model for implementing a 
partial model of a decision process of characteristics is 
attached and a result is substituted for a parameter of a 
5 linear model. The mechanistic model for altering a non- 
linear parameter value is established in form of a partial 
model independent from a main linear model, and there is a 
need to provide independently the associated mathematical 
model. With respect to the non-linear parameter, it is 

10 handled in form of the linear parameter in the linear model 
and it is set in a state that the linear logic is 
applicable. With respect to the non- linear model, 
necessary accuracy and reproducibility are varied in 
accordance with purpose of the use. Thus, there is a need 

15 to use a plurality of models in different ways even for the 
same non-linearity. Therefore, it is necessary for the 
me.chanistic model that a recombination of a model is 
permitted in a manner of nesting system. 
(2) Switch element 

20 With respect to implementation of a model of a 

structural non-linearity, it is possible to implement the 
model of the structural non-linearity in such a manner that 
an element, which serves as a switch for performing a 
connection and disconnection operation of the state 

25 quantity, is incorporated into the model. This switch 

element directly operates the state quantity of the model. 
Thus, it is possible to consider that the switch element is 
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one of parameters which determine a nature of the model. 
Of course, the switch element is dealt with as a variable 
which is the same as a co-efficient having values of "0" or 
"1" in the linear mathematical model. When the value is 
5 given by "1", the switch element serves to connect the 
state quantity, and when the value is given by "0", the 
switch element serves to disconnect the state quantity. 

This switch element needs a mechanistic model for 
a condition decision (formula) in which the predicted next 

10 state quantity is used to operate the switch. The decision 
is performed in such a manner that an operating quantity of 
"0" or "1" is generated from a predicted observation 
quantity predicted with the state of the model in the 
subsequent sampling time, and the switch element is 

15 operated (substituted) with this operation quantity in a 
manner that the internal connection is recombined before 
execution of the model through the operation of the switch 
element . Thus , the system and the associated mathematical 
model can be dealt with as the linearity even if the model 

20 is of the structural non-linearity. 

(3) Model including non-linearity 
From the above, a technique of implementing a 
model of non-linearity is classified as follows. One is 
that a mechanistic model is incorporated into a parameter 

25 included in a linear model so that the linear model is 

changed to a model which performs a non- linear behavior. 
In the transient characteristic of this model, the 
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characteristic value of the system has a continuity owing 
to a non- linear parameter which continuously changes, so 
the response is varied continuously. The other is that the 
switch element is included in the model of the structural 
non-linearity by a method in which a mechanistic model for 
predicting the subsequent state is incorporated to provide 
a non- linear structure for the model. In the transient 
characteristic of this model, the characteristic value of 
the system has discontinuity owing to connection and 
disconnection of the model structure by the switch element, 
so the response involves jumping phenomenon. 

From the above, the functional model is dealt with 
as the liner model regardless of linearity and non- 
linearity, and components to manufactured goods and up to 
the internal physical phenomenon can be modeled with 
greater accuracy and reproducibility. Further, as the 
mechanistic model is independent of a model to which a 
linear theory is applied, it is possible to incorporate 
thereinto many models utilized in the conventional product 
development, such as a statistical model, a recurrence 
model, an empirical formula, a fuzzy model and a neural 
network. In this manner, a non- linear model makes it 
possible to implement a model in a wider use, with which a 
non -linear theory, an engineering theory, a model, a stored 
technology, etc. are applied to such as structure, and 
mechanism common to products and components, physical 
phenomena and obstacle phenomena. Of course, it is 
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possible to reproduce all of those contents by an execution 
with a computer. 

(4) Implementation of nonlinearity of parameters 
It is known that a loss characteristic C 0 due to 
5 wind to a mobile car in running has a running resistance 
force f v proportional to the square of the velocity v d . 
That is, this is expressed by the equation f v =C 0 v d 2 . 
When this equation is expressed by the product of the 
velocity v d , it is possible to divide the equation into a 

10 linear equation of f v =Cv d and a non- linear loss 

characteristic C v =C 0 v d which is proportional to the 
velocity v d . In this relation, the equation is an 
equation representative of the entire system in form of 
the linearity, and the later is an equation representative 

15 of the non-linear loss characteristic. Thus, both the 

equations are associated with one another by an observation 
quantity v i =v d wherein v± is observed quantity of the 
velocity v d , and a coupling condition equation C = C V 
wherein C v is substituted for C . Figure 8 shows above- 

20 mentioned relation expressed by functional and mechanistic 
models . 

Fig. 1 shows a functional model in which input and 
output state quantity of running velocities v 1 • v 2 and 
driving resistance forces f i ' f 2 is provided with a loss 
2 5 characteristic C, and a mechanistic model in which an 

observation quantity v ± is put into a box to generate a 
non- linear loss characteristic C v . The mechanistic model 
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is expressed by a graph and a model. Part (a) of Fig. 1 is 
a graph in which the non- linear loss characteristic C v 
proportional to the absolute value of the observation 
quantity v ± of the velocity is expressed where abscissa 
(X axis) denotes the velocity v 1# and the ordinate (Y 
axis) denotes the non-linear loss characteristic C v . Part 
(b) of Fig. 1 is an example of the same content as the part 
(a) of Fig. 1 wherein the non-linear loss characteristic 
C v is expressed in form of a model with symbols shown in 
Figs. 78 and 79. Both the models are completely the same 
as the fundamental function element of the linear phase 
difference loss characteristic, except the mechanistic 
model of the non-linear loss characteristic C v . 
Incidentally, here, it is silent as to the estimated 
observation quantity is not mentioned, and its detail will 
be described later. 

- - The coupling condition of incorporating the 

mechanistic model into the functional model shown in Fig.l 
is given by the following mathematical model. 
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In eq. (1), the first line denotes a government 
equation of the functional model, the second line denotes a 
coupling condition equation for the observation quantity, 
and the third line denotes a coupling condition equation in 
which a loss characteristic is substituted. 

The mechanistic model incorporated into the 
functional model is given by the following mathematical 
model . 

C v = C 0 |v t . | ( 2 } 

It is possible to incorporate into the government 
equation of the first line of the eq. (1) the mathematical 
model of the mechanistic model, which determines the non- 
linear parameter of eq. (2), through the coupling condition 
equations on the second and the third lines of the eq. (1). 
(5) Modeling by the non-linear theory 
Let us consider further implementation of a model 
of the above-mentioned non-linear parameter. According to 
the example shown in Fig. 1, the non-linear loss 
characteristic C v proportional to the velocity v ± in eq. 
(12) is provided as the mechanistic model, and the equation 
is divided utilizing the fact that the loss of the running 
wind is proportional to the square of the velocity. 
However, such an example is rear on the actual non- linear 
model. Therefore, there is desired such a usually used 
general method that the non-linear model is transformed 



45 



directly to the mechanism model. As this idea, the non- 
linear model is divided or differentiated by the state 
quantity applied to the non-linear parameter before 
execution of the computation. Thus, it is possible to 
implement a device modeling of the mechanistic model . By 
way of example, representation of the model of Fig. 1 
according to this method offers the mechanistic model shown 
in Fig. 2. Incidentally, with respect to the dynamic 
frictional force F T shown in Fig. 2, it will be described 
later. 



Fig. 2 in form of a frame, part (a) of Fig. 2 shows a 
mathematical expression representation and part (b) of Fig. 
2 shows a model representation, wherein an estimated - 
observation quantity v ± _ at the subsequent sampling time 
of the velocity v d is applied to both the models. When 
-fv = C 0 V 1 _ 2 , which is an equation of force of a running 
wind resistance of the part (a) of Fig. 2, is generalized 
and expressed with a function f v = f n c (Vi_) , force f v of 
the running wind resistance of part (a) and part (b) of Fig. 
2 is expressed by the following mathematical model. 



With respect to the mechanistic models shown in 
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In eq. (3) 



the upper side is the force f v of the 
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running wind resistance expressed by the function wherein 
the equational expression in the parentheses denotes the 
mechanistic model, and the lower side is the equation of 
the estimated observation quantity v ± _ of the velocity V d . 
The estimated observation quantity v ± _ is an approximate 
value by a discrete system for estimating a velocity v d(k+1) 
at the subsequent sampling time using the velocity 
difference between the velocity v d(k) at the present 
sampling time and the velocity v d(k _ 1} at the previous 
sampling time. The absolute value of the estimated 
velocity Vi_ is the estimated observation quantity Vi_ of 
the input velocity v d of the non-linear loss 
characteristic C. Allowance of the application of the 
transformation of eq. (3) is limited to characteristics, 
factors and side loads. The method is the same also in the 
case that plural state quantities are inputted to the 
mechanistic model shown in Fig. 2. 

The mechanistic model shown in Fig. 2 can be 
represented by the following mathematical model. 
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In eq. (4), the upper side denotes the government 
equation, and the lower side denotes the coupling condition 
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equation between the functional model and the mechanistic 
model. The mechanistic model of the function represented 
by the upper portion of eq. (3) may be expressed by the 
following mathematical model. 
^ ^ 2 1 

(6) Modeling of structural non-linearity 
Next, the model of a tire kinetic frictional force 
F T appearing at the lower side of Fig. 2 is a model of a 

10 structural non-linearity in which a condition decision Swr 
predicts stop and run of a car from the estimated 
observation quantity v ± _ to operate a switch element Swr 
so that a side load of a kinetic frictional force F T0 "is 
controlled. In this structural non-linearity, when the 

15 velocity Vi_ is 0, the condition decision S WT is 1 so that 
th-e switch element Swt turns off, and thus the side load 
F T0 is disconnected. On the other hand, when the velocity 
v ± _ is 0, the side load F T0 is connected. The result is 
added to the force f v of the running wind resistance in 

20 form of the tire kinetic frictional force F T . 

The condition decision equation and the tire 
kinetic frictional force are expressed as follows. 

/O; _ = 0) theniS^ = 1) else (5^ = 0)' 

F = S F I ( 6 ) 

25 In eq. (6), the upper side denotes the condition 
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decision equation, and the lower side denotes the tire 
kinetic frictional force which is the same as one 
incorporated into the equation appearing at the upper side 
of eq. (4) . 

5 Here, there will be supplemented a method of 

condition decision taking a decision of the zero velocity 
in eq. (13) by way of example. In the condition decision, 
it often happens that a zero value or and a designation 
value of the state quantity are provided for a decision. 

10 It is difficult for the actual simulation including a 

calculation error to determine mathematically this decision. 
In some cases, a failure of determination causes an 
oscillation of the state quantity in synchronism with a 
sampling period. For example, in eq. (6), even if a 

15 velocity within a certain range is set up to be zero so 

that the calculation error can be avoided, there will occur 
an oscillation phenomenon owing to an alteration of the 
sampling period, or increment or decrement of the 
acceleration. As one of methods of avoiding this problem, 

20 there is a method of predicting and deciding a stability of 
a velocity to zero on the basis of the velocity v d(k) at 
the present time (k), the velocity difference A v d from the 
previous sampling time (k-1), and the velocity v d(k+1) at 
the subsequent sampling time (k+1). This condition 

25 decision is expressed as follows. 



49 



if ~^h|Av|) then (S^ =1) else (S^ = 0) 

In eq. (7), the upper side denotes the velocity 
difference from the previous sampling time, and the lower 
side denotes the condition decision equation, where v d(k) 
5 is a velocity at the present time; V TH , a criterion value; 
v d(k-i) / a velocity at the previous sampling time; and A v , 
a velocity difference from the previous sampling time. 
With respect to Swt appearing at the lower side, when a 
condition is satisfied in a decision condition of the 

10 velocity 0, S WT = 1 , or else S m =0 . Incidentally, it is 

effective that the criterion value for the zero velocity is 
given by V TH = 0 . The algorithm for a decision decides 
whether the velocity difference A v is associated with a 
state that the criterion value V TH is intervened between 

15 the velocity v d(k) at the present time and the velocity V 
d(k+i> at "t ne subsequent sampling time. That is, it is 
decided whether the velocity v d(10 at the present time is 
within a range wherein the velocity difference A v is 
established as the criterion width. Consequently, 

20 according to the condition decision of operating a switch 
element while a velocity stop and a convergence state to a 
target value are decided, since those are converged while 
the decision width is sequentially reduced, it is possible 
to obtain the convergent state less in the residual between 

25 the criterion value and the decided state quantity. And 

also it is possible to converge an oscillation occurring in 
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synchronism with a sampling period taking the decided value 
as a boundary. 

For the condition decision, there is used usually 
a method in which some latitude is allowed for the 
criterion value to stabilize the decision result. In the 
event that some latitude is allowed for the criterion in 
the condition decision appearing at the lower side of eq. 
(7), the following equation is given. 

if (S^ = 0) n (\v d(k) -V SH \< | Av|) then (S^ = 1) 1 
else if (5 W =1) f](\v d(k) -V SL \< |Av|) then (S^ = 0) J 

(8) 

In eq. (8), V SH denotes the upper limit of the 
criterion, and V SL denotes the lower limit of the 
criterion. The remaining is the same as eq. (7). With 
respect to the decision result of the eq. (8), in the event 
that the velocity is smaller than the upper limit V SH of 
the criterion (velocity up), when the velocity reaches the 
upper limit value V SH or larger, the decision result is 
given by S WT = 1 . And, in the event that the velocity is 
larger than the lower limit V SL of the criterion (velocity 
down), when the velocity reaches the lower limit value V SL 
or smaller, the decision result is given by S WT =0. 
(7) Inappropriate model representation 
Next , there will be given an example of device 
modeling in which a non- linear parameter is not appropriate. 
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Fig. 3 shows an example of device modeling wherein 
the non-linear loss characteristic C shown in Fig. 1 is 
given directly as a non-linearity with use of the state 
quantity in a functional model. The government equation of 
5 this model is given by the following expression. 

fl =CvA^) + f 2 

The functional model in Fig. 3 becomes eq. (9) 
which directly operates the state quantity, and becomes an 
unusual mathematical model representative of a non- 
10 linearity inherent to a running resistance. Thus, it is 
difficult to represent such a mathematical model by a 
linear government equation. Further, regarding products 
and components having many types of non-linearity as their 
inherent functions, it is almost impossible to solve their 
15 coupling physical phenomena directly with the mathematical 
model. This is considered as a limit of the conventional 
non- linear model. 

For example, comparing eqs.(l) and (2) with eq. 
(9), eqs.(l) and (2) are divided into the expression parts 
20 of the non-linearity and of the linearity, even if the same 
model is concerned. Thus, it is possible to solve the 
equation through a substitution of each state quantity in 
plural linearity expressions. On the other hand, in case 
of the eq. (9), the expression part of the non-linearity 
25 and the expression part of the linearity are mixed in one 
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equation. Thus, it is difficult to solve the equation 
through a substitution of each state quantity of plural 
non-linearity expressions. 

In this manner, for modeling of the device with 
the non-linearity, there are several ways of methods of 
expression even if function and behavior are same. However, 
it would be understood that direct insertion of non-linear 
model into linear model limits device modeling. In order 
to avoid this limitation, it is important that the 
functional model and the' government equation are modeled in 
form of the linearity and the non-linearity to be inserted 
into this functional model is modeled separately as the 
independent mechanistic model. The device modeling, which 
fails to satisfy this condition, is associated with the 
following problems. 

CD It is difficult to solve mathematically 
products goods and components in which plural of non- 
linearities are involved. 

CD A linear theory cannot be applied to the 
government equation. 

CD A model of a non- linear parameter cannot be 
recombined in a manner of nesting system. 

A functional model and a model with a non- 
linear parameter cannot be standardized. 

From the above, the model shown in the parts (a) 
and (b) of Fig. 10 is an example of a device modeling, 
which is inappropriate as the functional model. 
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However, with respect to the representation shown 
in the parts (a) and (b) of Fig. 3, the following uses are 
applicable with a consideration for processing the 
representation on a software basis. 
5 CD The uses of a flow diagram and an explanatory 

view for simplifying a complicated system to visually 
understanding of the whole. 

(D Expression means for simplifying a functional 
model used repeatedly in a development field. 
10 (3) Symbols for symbolizing functional models and 

parameters including a non-linearity. 

(8) Running resistance model of a car 

The functional model shown in Fig. 2 is an example 
in which a structural non-linear model of the dynamic 
15 frictional force F T is added to the non-linear loss 
characteristic C v shown in Fig. 1 so that a running 
resistance is improved in accuracy. The non-linear loss 
characteristic C v and the dynamic frictional force F T of 
the tire shown in Fig. 2 are running loads which are 
20 important in deciding running performance and fuel 

efficiency of a mobile car. More in detail, this may be 
expressed by the following equation. 

f T =^pC d Av d 2 + v m W m cos(6)+W m sin(0) 

25 (1 0) 

where p denotes air density; C d coefficient of 
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air resistance; A frontal projected area; fi m coefficient 
of rolling resistance; W m weight of vehicle; and 9 an 
angle of inclination. In this equation, the first term 
denotes a force f v by the running wind, the second term 
5 denotes a dynamic frictional force F T of a tire, and the 
third term denotes a force of an up-and-down slope by a 
weight. When the first and second terms of the eq. (10) 
are rearranged with the absolute value of the velocity v d 
after the example of the eq. (3) to make up the whole into 

10 one mechanistic model in 1 form of equivalent non-linear loss 
characteristic C T . Functional and mechanistic models, 
wherein mass W m / g of a mobile car is added to this 
mechanistic model, is shown in Fig. 4 where v d denotes 
velocity of the car; f d driving force; f T force of 

15 running resistance; and x m internal state quantity. 

Fig. 4 shows a mechanitic model of equivalent non- 
linear loss characteristic C T in which a hierarchical 
structure developed below is expressed by a relation of 
inclusion. This model incorporates into the lower portion 

20 thereof a mechanistic model for reproducing a potential 
quantity of the running wind load f v and the kinetic 
frictional force F T so as to form the mechanistic model 
in which the running wind load f v and the dynamic 
frictional force F T are added and then combined to the 

25 upper portion of the model and the equivalent non- linear 
loss characteristic C T , which is transformed after the 
example of eq. (10), is substituted for C. It is 
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understood that the mechanistic model of Fig. 4 is 
developed vertically and horizontally into systematic 
device modeling. 

The government equation and the coupling condition 
equation of the mechanistic model of Fig. 4 are as follows. 



w m 
g 

0 



-C 1 -W m sm{9) 
-10 0 



(11) 



c = c v 

The mechanistic model of Fig. 4 can be expressed 
by the following mathematical model. 



J4smp , 
„ X - 



dt 



/v 



\pC d Av it , 2 



C T =(f v + F T ) 1 ± 1 



(12) 



In eq. (12), the first line denotes an estimated 
velocity v L ; the second line a coupling condition equation 
for distributing the entered estimated velocity v ± to the 
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subordinate mechanistic model; the third line an equation 
for reproducing a force f v of a running wind resistance; 
the fourth line a dynamic frictional force F T of a tire; 
and the fifth line an equivalent non-linear loss 
5 characteristic C T . With respect to variables of eq. (12), 
v z denotes an observation quantity of a velocity v d ; Swt 
a condition decision of the upper side of eq. (6); v iv a 
velocity entered to a mechanistic model of a running wind 
resistance developed on a subordinate basis; and v iT a 

10 velocity entered to a mechanistic model of a tire 

frictional resistance developed on a subordinate basis . 

As shown in Fig. 4, the mechanistic model can be 
unified to a single mechanistic model, including the 
structural non-linearity including a switch element. "This 

15 shows that as far as the non-linearity of the same function 
is concerned, it is regarded as one system, and it is 
possible to systematize the system through a vertical 
development in its entirety. 

Hereinafter, there will be described embodiments 

20 of the present invention. 

According to the embodiments of the present 
invention, described hereinafter, each of the embodiments 
comprises a computer and a program for causing the computer 
to operate as a non-linear characteristic reproducing 

25 apparatus according to one embodiment of the present 

invention. First, the computer will be explained, and then 
the function of the non- linear characteristic reproducing 
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apparatus according to one embodiment of the present 
invention, which is implemented by the computer, will be 
explained . 

Fig. 5 is a perspective view of a computer in 
5 which a non-linear characteristic reproducing apparatus 
according to one embodiment of the present invention is 
implemented . 

A computer 10 comprises: a main frame 11 
incorporating thereinto a CPU, a main memory, a hard disk, 
10 a communication board, etc.; a display unit 12 for 

displaying images and a character string on a display 
screen 12a in accordance with the main frame 11; a key 
board 13 for inputting a user's instruction to the computer 
10; and a mouse 14 for designating an optional position on 
15 the display screen 12a to enter an instruction according to 
an icon or the like displayed on the position at the time 
of. the designation. 

The main frame 11 of the computer 10 is provided 
with an FD loading slot 11a and a CDROM loading slot lib 
20 for loading a floppy disk and a CDROM (not illustrated in 
Fig. 5, see Fig. 6), respectively. Inside the main frame 
11, there is incorporated a floppy disk drive and a CDROM 
drive for driving and accessing the floppy disk and the 
CDROM loaded through the FD loading slot 11a and a CDROM 
25 loading slot lib, respectively. 

Here, the CDROM loaded onto the computer 10 stores 
the non-linear characteristic reproducing program referred 
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to the present invention. The CDROM is loaded into the 
main frame 11 through the CDROM loading slot lib so that 
the non- linear characteristic reproducing program stored in 
the CDROM is installed into a hard disk of the computer 10 
5 by the CDROM drive. When the non -linear characteristic 
reproducing program installed into the hard disk of the 
computer 10 is started, the computer 10 operates as one 
embodiment of the non-linear characteristic reproducing 
apparatus according to the present invention. 

10 Therefore, the CDROM storing the non-linear 

characteristic reproducing program corresponds to one 
embodiment of the non -linear characteristic reproducing 
program storage medium according to the present invention. 

The non- linear characteristic reproducing program 

15 stored in the CDROM is installed into a hard disk of the 

computer 10 in the manner as mentioned above. In this case, 
the hard disk, into which the non-linear characteristic 
reproducing program is installed, also corresponds to one 
embodiment of the non-linear characteristic reproducing 

20 program storage medium according to the present invention. 
Further, in the event that the non -linear 
characteristic reproducing program is down loaded onto a 
floppy disk and the like, the floppy disk and the like, 
which store the non-linear characteristic reproducing 

25 program down loaded, also correspond to one embodiment of 
the non- linear characteristic reproducing program storage 
medium according to the present invention. 



59 



Fig. 6 is a schematic diagram showing a hardware 
structure of the computer shown in Fig. 5. 

Here, the computer comprises a central processing 
unit (CPU) 111, a RAM 112, a hard- disk controller 113, a 
5 floppy disk drive 114, a CDROM drive 115, a mouse 

controller 116, a keyboard controller 117 and a display 
controller 118. They are mutually connected via a bus 110. 

The floppy disk drive 114 and the CDROM drive 115, 
onto which the floppy disk 51 and the CDROM 520 are loaded, 
10 respectively, as mentioned above referring to Fig. 5, 

access the floppy disk 51 and the CDROM 52, respectively. 

Here, there are further shown a hard disk 53 to be 
accessed by the hard-disk controller 113, a mouse 14 to be 
controlled by the mouse controller 116, a keyboard 13 to be 
15 controlled by the keyboard controller 117, and a display 
unit 12 to be controlled by the display controller 118. 
,- - Each of Figs. 7 to 11 is a typical illustration 

showing a structure of a non-linear characteristic 
reproducing program stored in a storage medium. 
20 In Figs. 7 to 11, storage media 200, 300, 400, 500 

and 600, are representative of the above-mentioned CDROM, 
hard disk, floppy disk and the like. 

Fig. 7 shows an embodiment of a first non-linear 
characteristic reproducing program storage medium. The 
25 storage medium 200 stores a non-linear characteristic 
reproducing program 210 having a state quantity 
transforming section 220 and a non- linear characteristic 
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reproducing section 230. The non-linear characteristic 
reproducing program 210 causes the computer to serve as a 
non-linear characteristic reproducing apparatus which 
receives a predetermined first state quantity and outputs a 
5 second state quantity by applying the non-linear 

characteristic reproducing process to the first state 
quantity. 

Here, the a state quantity transforming section 
220 performs a linear formation from the first state 

10 quantity to the second state quantity for each sampling 

time in accordance with a transformation parameter set up. 
The non-linear characteristic reproducing section 230 
determines a transformation parameter for transformation of 
at least one of the first state quantity and the second 

15 state quantity, or a state quantity derived from at least 
one of the first state quantity and the second state 
quantity at the subsequent sampling time in accordance with 
the estimated observation quantity at the subsequent 
sampling time, and sets the transformation parameter thus 

20 determined to the state quantity transforming section 220. 
In some case, the non-linear characteristic reproducing 
section 230 may receive the estimated observation quantity 
and one or more variables as well and determine a 
transformation parameter in accordance with the entered 

25 estimated observation quantity and one or more variables. 

Specifically, in the non-linear characteristic 
reproducing section 230, a standardized estimation value. 
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in which an estimation value of the second state quantity 
at the subsequent sampling time is standardized with the 
above-mentioned estimated observation quantity, is 
determined in form of the transformation parameter. More 
5 specifically, in the non-linear characteristic reproducing 
section 230, an estimation value of the second state 
quantity at the subsequent sampling time is subjected to a 
division or a differential with the absolute value of the 
above-mentioned estimated observation quantity, so that the 
10 estimated observation quantity can be determined in form of 
the transformation parameter. 

This non- linear characteristic reproducing section 
230 is, as will be described later, preferably applied to 
reproductions of, for example, a characteristic of a non- 
15 linear spring, a characteristic of an air spring, a 

characteristic of a link mechanism, a characteristic of an 
object which moves involving friction, and a characteristic 
of a mechanism for variable moment of inertia in which 
moment of inertia is varied in an angular velocity. 
20 Fig. 8 shows an embodiment of a second non -linear 

characteristic reproducing program storage medium. A 
storage medium 300 stores therein a non-linear 
characteristic reproducing program 310 having a linear 
model section 320 and a non-linear model section 330. The 
25 non-linear model section 330 constituting the non-linear 
characteristic reproducing program 310 has non- linear 
transforming sections 331a, 331b ••• 331n. The number of 
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non-linear transforming sections is according to the way. 
In some case, it happens that the non-linear characteristic 
reproducing program has only one non-linear transforming 
section. Fig. 8 shows a general case of the use of a 
5 plurality of non- linear transforming sections. 

The non-linear characteristic reproducing program 
310 causes the computer to serve as a non -linear 
characteristic reproducing apparatus for reproducing 
characteristics of systems including a non-linear system. 

10 The linear model section 320 reproduces characteristics of 
a linear system. The non- linear model section 330 obtains 
from the linear model section 320 an estimated observation 
quantity at the subsequent sampling time with respect to at 
least one of a first state quantity and a second state 

15 quantity, which are mutually in a relation of a non -linear 
transformation, or a state quantity derived from said at 
least one of the first state quantity and the second state 
quantity, and determines a transformation parameter used in 
a linear transformation at the subsequent sampling time 

20 between the first state quantity and the second state 

quantity. An operation for determining the transformation 
parameter at the subsequent sampling time by the non-linear 
model section 330, and a linear operation including a 
linear transformation between the first state quantity and 

25 the second state quantity using the transformation 

parameter at the subsequent sampling time obtained by the 
non- linear model section 330 are alternately repeated. 



63 



The non-linear transforming sections 331a, 331b ••• 
331n, which constitute the non-linear model section 330, 
are portions dealing with a first state quantity and a 
second state quantity which are mutually in a relation of a 
5 non- linear transformation and are different in sort of a 

state quantity of at least one of the first state quantity 
and the second state quantity. Each of the non- linear 
transforming sections 331a, 331b "• 331n obtains an 
estimated observation quantity at the subsequent sampling 

10 time with respect to at least one of the first state 

quantity and the second state quantity, which are dealt 
with by the associated non-linear transforming section, or 
a state quantity derived from said at least one of the 
first state quantity and the second state quantity, and 

15 determines a transformation parameter used in a linear 

transformation at the subsequent sampling time between the 
first state quantity and the second state quantity. 

The non- linear characteristic reproducing program 
310 is, as will be described later, preferably used in, for 

20 example, reproducing of characteristics of Geneva mechanism 
and reproducing of characteristics of a liquid residue 
warning lamp wherein a lamp is connected in series with a 
thermistor. 

Fig. 9 shows one embodiment of a third non-linear 
25 characteristic reproducing program storage medium. A 
storage medium 400 stores therein a non -linear 
characteristic reproducing program 410 having a logical 
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decision section 420 and. a state quantity selection section 
430. The non-linear characteristic reproducing program 
causes the computer to serve as a non-linear characteristic 
reproducing apparatus for reproducing characteristics of 
5 systems including a non- linear system. 

The logical decision section 420 receives one or 
more variables and determines a logical value at the 
subsequent sampling time, which is selected from among a 
plurality of discrete values, in accordance with the one or 
10 more variables entered. 1 

The state quantity selection section 430 receives 
a predetermined input state quantity and outputs an output 
state quantity in which at the subsequent sampling time, a 
relation between the input state quantity and the output 
15 state quantity is switched into a relation according to the 
logical value at the subsequent sampling time which is 
determined by the logical decision section 420. 

In some case, the state quantity selection section 
430 outputs at the subsequent sampling time an output state 
20 quantity in which a connection relation between the input 
state quantity and the output state quantity is switched 
into a connection relation according to the logical value 
at the subsequent sampling time which is determined by the 
logical decision section 420. Or alternatively, in some 
25 case, the state quantity selection section 430 outputs an 
integrated input state quantity in form of an output state 
quantity. In this case, in the event that the logical 
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value at the subsequent sampling time, which is determined 
by the logical decision section 420, is a predetermined 
logical value, the state quantity selection section 430 
outputs at the subsequent sampling time an output state 
5 quantity which is changed over to the initial value. 

Further, in some case, the logical decision 
section 420 determines a logical value at the subsequent 
sampling time selected among from a plurality of discrete 
values in accordance with an estimated observation quantity 

10 at the subsequent sampling time for a plurality of input 

state quantities. And the state quantity selection section 
430 receives the plurality of input state quantities and 
outputs at the subsequent sampling time as an output state 
quantity an input state quantity selected according to the 

15 logical value at the subsequent sampling time determined by 
the logical decision section 420. 

- - This non- linear characteristic reproducing program 

410 is, as will be described later, preferably applied to 
reproductions of, for example, a characteristic of a 

20 mechanism including two members having a relative movement 
possible state and a relative movement impossible state, a 
characteristic of a clutch mechanism including two members 
having a slipping state of relatively slipping and a 
connecting state of mutually connecting, a characteristic 

25 of a brake mechanism for applying a braking energy to a 
drive shaft, and a characteristic of an automatic-reset 
mechanism in which when a spring is effected on a movable 
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member so that an application of the external force to the 
movable member is released, the movable member is 
automatically reset to the initial state by the effect of 
the spring, said automatic-reset mechanism being provided 
5 with a stopper for restricting a movable range of the 
movable member. 

Fig. 10 shows one embodiment of a fourth non- 
linear characteristic reproducing program storage medium. 
A storage medium 500 stores therein a non-linear 

10 characteristic reproducing program 510 having a state 
variation estimation section 520, a state deviation 
detection section 530, a stable state decision section 540 
and a state quantity selection section 550. 

The non- linear characteristic reproducing program 

15 causes the computer to serve as a non-linear characteristic 
reproducing apparatus for receiving a predetermined first 
state quantity and outputting a second state quantity in 
which the second state quantity is changed over to a non- 
linearity of one according as the first state quantity is 

20 in a predetermined state. 

The state variation estimation section 520 
performs a prediction of a state quantity variation width 
of the first state quantity during a period from the 
present sampling time to the subsequent sampling time. 

25 The state deviation detection section 530 

determines a deviation between a value of the first state 
quantity at the present sampling time and a decision value 
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for determining whether the first state quantity is in a 
predetermined state. 

The stable state decision section 540 determines a 
logical value, which is selected among from a plurality of 
5 discrete values, predicting a variation of the non- 
linearity at the subsequent sampling time in accordance 
with size between the state quantity variation width 
predicted by the variation estimation section 520 and the 
deviation determined by the state deviation detection 
10 section 530. 

The state quantity selection section 550 outputs 
the second state quantity changed over in accordance with 
the logical value determined by the stable state decision 
section 540. 

15 Here, in some case, the state deviation detection 

section 530 determines a deviation between a value of the 
first state quantity at the present sampling time and a 
decision value for determining whether the first state 
quantity is in a predetermined state, the decision value 

20 being varied in accordance with positive and negative of 
the state quantity variation width predicted by the 
variation estimation section 520. 

Fig. 11 shows one embodiment of a fifth non- linear 
characteristic reproducing program storage medium. A 

25 storage medium 600 stores therein a non-linear 

characteristic reproducing program 610 having a linear mode 
section 620 and a non-linear model section 630. The non- 



68 



linear characteristic reproducing program 610 causes the 
computer to serve as a non- linear characteristic 
reproducing apparatus for reproducing characteristics of 
systems including a non- linear system. 
5 The linear mode section 620 reproduces 

characteristics of a linear system, and includes a state 
quantity transformation section 621 for linear -transforming 
a first state quantity to a second state quantity in 
accordance with the transformation parameter set up. 

10 The non-linear model section 630 obtains from the 

linear mode section 620 an estimated observation quantity 
at the subsequent sampling time of a predetermined 
observation state quantity to generate a transformation 
parameter in accordance with the estimated observation 

15 quantity, and sets the transformation parameter thus 

generated to the state quantity transformation section 621. 
The non- linear model section 630 comprises a slow change 
reproducing section 631 and a characteristic generation 
section 632. According to the present embodiment, the slow 

20 change reproducing section 631 comprises a stationary value 
set section 631a and a normalization response section 631b. 

The slow change reproducing section 631 obtains 
from the linear mode section 620 observation quantity or 
estimated observation quantity of a predetermined second 

25 observation state quantity, which is the same as the first 
observation state quantity or different from the first 
observation state quantity, and generates a slow change 
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state quantity reflecting characteristics of a first non- 
linear system representative of a relatively slow behavior 
change in accordance with the observation quantity or 
estimated observation quantity of the second observation 
state quantity. 

Of the slow change reproducing section 631, the 
stationary value set section 631a determines a stationary 
value of the slow change state quantity after infinite time 
elapsed in the event that the state of the linear mode 
section 620 is continued, in accordance with the 
observation quantity or estimated observation quantity of 
the second observation state quantity. And the 
normalization response section 631b generates the slow 
change state quantity reflecting characteristics of a first 
non-linear system, which is to be transferred to the 
characteristic generation section 632, in accordance with 
the stationary value of the slow change state quantity 
determined by the stationary value set section 631a, and 
the known normalized time variation characteristic of the 
slow change state quantity in the first non-linear system. 

The characteristic generation section 632 
generates a transformation parameter reflecting 
characteristics of a second non- linear system 
representative of relatively rapid behavior change in 
accordance with the estimated observation quantity of the 
first observation state quantity at the subsequent sampling 
time, which is obtained from the linear mode section 620, 
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and the slow change state quantity generated by the slow 
change reproducing section 631, and then sets the generated 
transformation parameter to the state quantity 
transformation section 621. 
5 The non- linear characteristic reproducing program 

610 is preferably used, for example, in reproductions of 
characteristics of a system having a device which varies in 
resistance value in accordance with temperature variation. 

Next, there will be explained various functions as 
10 a non- linear characteristic reproducing apparatus of the 
present invention, which is implemented in the computer 
shown in Fig. 5, when the non-linear characteristic 
reproducing programs shown in Figs. 7-11 are installed in 
the computer. 

15 For functional blocks, which will be explained 

hereinafter, the same names as functional blocks (software 
components) on a software shown in Figs. 7-11 are used. It 
is noted, however, that each of the functional blocks as 
the non-linear characteristic reproducing apparatus, which 

20 will be explained hereinafter, denotes a complex of a 
functional block on a software for implementing the 
associated function and a structural element on a hardware. 
On the other hand, each of the functional blocks on a 
software shown in Figs. 7-11 denotes a software (an 

25 application software), excepting a hardware, an operation 
system, etc., of elements for implementing the associated 
function . 
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Fig. 12 is a principle explanatory view of a first 
non-linear characteristic reproducing apparatus according 
to an embodiment of the present invention. Incidentally, 
Fig. 12 is a view useful for understanding the principle of 
5 the first non-linear characteristic reproducing apparatus 
of the present invention, and the first non-linear 
characteristic reproducing apparatus of the present 
invention is not restricted to Fig. 12 and the associated 
description . 

10 In the right side of Fig. 12, a state quantity 

transformation section R denotes characteristic values such 
as a coefficient of resistance, an electric capacity, mass, 
a coefficient of link-transf er , and input and output state 
quantities v and f denote state quantities such as voltage, 

15 velocity, force and current. The state quantity 

transformation section R transfers the input state quantity 
v -to the output state quantity f. The non-linear 
characteristic reproducing section appearing at the right 
side of Fig. 12 receives an estimated observation quantity 

20 v_ (an outline half -circle symbol on a colored background) 
at the subsequent sampling time of the input state quantity 
v applied to the state quantity transformation section, and 
in addition external state quantities v L and f i and an 
external signal S x . These input state quantities cause 

25 the non- linear characteristic reproducing section to 

generate an estimated value f_= (v_, S l9 f i) of the 

output state quantity at the subsequent sampling time. The 
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output state quantity is divided or differentiated with the 
estimated observation quantity v_ to generate a prediction 
value of characteristic R p of the state quantity 
transformation section at the sampling time. The 
5 characteristic R p is substituted for the state quantity 
transformation section R (an arrow of the dot line), and 
then a calculation for the whole system model into which 
the state quantity transformation section is incorporated 
is carried out. Repeatedly executing this operation for 

10 each sampling time makes it possible to reproduce non- 
linear physical phenomena and behavior of products and 
components in form of time recording data. 
(Explanation of examples) 
1. Flow storage element (Conical spring) 

15 (1) Functional model 

Fig. 13 shows an example of a mechanistic model in 
which with respect to a linear flow storage element 
(fundamental functional element), stiffness of a spring 
varying according to quantity of deformation is given with 

20 a graphical representation. Fig. 13 shows the linear 

mechanistic model where v x and v 2 denote velocity, f i and 
f 2 denote an external force, and x denotes an internal 
state quantity. v ± denotes a deformation velocity 
differential of the spring. The mechanistic model 

25 incorporated into Fig. 13 is a graph showing that stiffness 
K d of the spring is a function of a quantity of deformation 
d . This graph can represent a linear spring to a non- 
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linear spring. 6 0 in Fig. 13 denotes a quantity of 
deformation generating an offset load. 

A flow storage element characteristic can be 
expressed by the following government equation and coupling 
5 condition equation. 
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(13) 

/ 2 J 

K = K d 

The mechanistic model of stiffness of the spring 
shown in Fig. 13 can be expressed by the following 
mathematical model. 



d =jv d dt 



K d =fric(3 +3 0 )\ 



(14) 



In Fig. (14), the first line denotes a deformation 
6, and fnc (a) on the second line denotes a function 
15 determining stiffness K d of a non-linear spring. 

(2) Mechanistic model (non- linear spring model) 
As a typical example of a non-linear spring, there 
is a conical spring. Part (a) of Fig. 14 shows a structure 
of a conical spring wound at constant pitch. Part (b) of 
20 Fig. 14 shows a graph representative of a relation between 
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deformation 6 and load F. Regarding symbols in the figure, 
H 0 denotes a free length; Rj. an effective radius on a 
small side; R 2 an effective radius on a large side; d a 
diameter of a coil; and P a pitch of the coil. The spring 
5 is defined with winding number n of a free length and 
modulus G of transverse elasticity. 

According to the conical spring shown in the part 
(a) of Fig. 14, a torsion moment on the side associated 
with a larger diameter is large and it is more easily 

10 deformed, and thus a contact of the coils is sequentially 

initiated from the end larger in diameter. In the part (b) 
of Fig. 14, regarding the spring characteristic, a relation 
between deformation <5 and load F offers proportion during 
a range 0 to A before the contact of the coils, and the 

15 effective winding number is reduced during a range A to B 
after the contact of the coils, so that the load F to the 
deformation 6 is increased. And finally, on the point B 
the smaller diameter end contacts to offer a bottom- 
projection so as to form rigid body. 

20 With respect to this conical spring, a relational 

expression is referred from the well known reference. 
Between the load F and the deformation <5 , the following 
relation is applied. 

< Spring load and deformation before contact 

25 (linear area) > 
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d = 



16n 
Gd 2 



(15) 



< Spring load and deformation after contact (non- 
linear area) > 



R = R 2 -(R 2 -R 1 )- 



d'=dJl 
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Equation (16) represents a non- linear area after 
contact. In eq. (16), the first line denotes an 
approximate expression of an average radius R when the n' - 
th coil from R 2 side contacts, the second line denotes a 
central distance d' in the vertical direction between the 
contacted coils, and the third and fourth lines denotes a 
relation between the load F and the deformation (5 . From 
equations of the conical spring of eqs.(15) and (16), a 
device modeling is performed on the mechanistic model to be 
incorporated into Fig. 13. 

The respective expressions of eq. (16) are 
rearranged and the linear expression of eq. (15) is given 
in form of a non-linear expression with a load ratio O . In 
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this case, the following equations are provided. 



5 




(17) 



(18) 



(19) 



Equation (17) represents load-deformation 
characteristics of the conical spring. With respect to 
stiffness K of the spring in the linear area, R of eq. 
(18) is R 2 , and the load ratio of eq. (18) is given by <E> 

15 = 1 . In the non-linear area, effective radiuses R are 
sequentially determined for each contacted coil from the 
first line of eq. (16), and the effective radiuses R thus 
determined is substituted for R of eqs.(18) and (19). A 
relation between the load F to the deformation d , which is 

20 determined in accordance with the above-mentioned method, 
is previously computed, and stiffness K d of the spring, 
which is determined through deformation <5 retrieval by 
integration of estimated velocity v<3_ , is substituted for 
K of the functional model. The functional and mechanistic 

25 models representative of this relation is shown in Fig. 15. 
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Fig. 15 shows one in which the mechanistic model 
of the conical spring is substituted for the functional 
model of Fig. 13. The functional on the second line of eq. 
(14) is replaced by a function to stiffness K of the spring 
5 from the deformation 6 . 

2 . Example of device modeling of an air spring 
By way of example of a fluid system, there is 
performed a device modeling of a pneumatic cylinder 
simulating a compression and expansion process of a car 
10 engine. Incidentally, modeling of the internal energy 

(temperature energy) caused by compression and expansion is 
omitted. 

In Fig. 16, p 0 denotes an atmospheric pressure; 
S a cross-sectional area of a cylinder; M mass of a 

15 piston; C a coefficient of viscosity resistance; K 

stiffness of an air spring; D a coefficient of internal 
damping; f a piston external force; v a piston velocity; 
P a cylinder internal pressure; and V a cylinder volume. 
In Fig. 16, the external force f is effected on a piston, 

20 so that the piston velocity v and the associated cylinder 
volume change A V p occur whereby the cylinder internal 
pressure p generates. The cylinder internal pressure p 
effects as the reaction force in form of stiffness K of an 
air spring. The rise of the internal pressure involves the 

25 temperature rise inside the cylinder, so that the thermally 
expanded volume increase component becomes a part of the 
internal pressure rise component . The raised temperature 
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is radiated in accordance with the coefficient of internal 
damping D so that the thermally expanded volume increase 
component is decreased. Heat insulation and heat radiation 
of the cylinder is effected in accordance with the 
5 coefficient of internal damping D. Incidentally, the 

initial state of the cylinder volume is set to a cylinder 
volume V 0 in which the internal pressure p is equal to 
the atmospheric pressure p 0 . 

(1) Functional model 

10 Compression and expansion process of the piston, 

which summarizes this relation, is represented by a 
functional model shown in Fig. 17. Heat insulation and heat 
radiation of the cylinder is effected in accordance with 
the coefficient of internal damping D. According to the 

15 functional model shown in Fig. 17, a mechanistic model of 
an air spring is incorporated into a flow storage element 
f fundamental functional element). A mechanistic model of 
temperature rise will be omitted. 

Fig. 17 is one for determining a cylinder volume 

20 (an estimated accumulation observation quantity) where P 

denotes an observation quantity of an internal pressure of 
a cylinder; V P _ a double integration mark; and x aa a 
movement distance of a piston. 

In Fig. 17, the government equation, the coupling 

25 condition and the estimated observation quantity are given 
by the following equation in order. 
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K = K, 



(2 0) 



In eq. (20), on the upper side, first and second 
lines denote a state equation; the third line a state 
equation of a storage observation quantity where x aa 
denotes a movement distance (a storage potential quantity) 
of a piston; the fourth line an input and output equation; 
and the fifth line an equation of an observation quantity 
of pressure in a cylinder. The lower side denotes a 
coupling condition equation to be connected with a 
mechanistic model where K s denotes stiffness of an air 
spring in the mechanistic model, and Vi_ denotes volume of 
a movement distance of a piston. 

(2) Mechanistic model (air spring model) 

Modeling of the mechanistic model shown in Fig. 17 
is as follows. 

First, the cylinder volume V determined by a 
movement quantity of the piston shown in Fig. 17 is 
expressed by the following equation. 
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V p_ = ^ SX aa -j^Xc^smp 

V=V V -(Sx aa +AV p _) 



(2 1) 



In eq. (21), the upper side denotes an estimated 
storage observation quantity of a volume change, and the 
lower side denotes an estimated observation quantity of 
volume of the cylinder. Incidentally, t smp denotes a 
sampling period. 

Next, a relation between the pressure p in a state 
of compression and the volume V is expressed in form of the 
polytropic change by the following equation where n is an 
index . 

P = P °{^Fj (22) 

In eq. (22), stiffness K s of an air spring is 
varied as the index n determining the polytropic change in 
accordance with the following selection of the index n. 



n= 0 
1 



isobaric change 
isothermal change 

adiabatic change ( K, is ratio of specific heats) 
isovolume change 



The non- linear parameter, in which eq. (22) is 
transformed, is expressed by the following equation. 
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Ap c =x\- 



(2 3) 



In eq. (23), the upper side denotes a linear 
equation representative of a relation between a cylinder 
5 pressure variation Ap and a volume change AV p _with an 
air spring stiffness K s , and the lower side denotes a 
non-linear equation of an air spring stiffness K s in the 
volume V from eq. (22). Both the equations causes the 
model to offer a relation in which the linearity and the 

10 non-linearity are separated and are independent one another. 

The air spring stiffness K s can be determined by 
means of differentiating eq. (23) with the volume change. 
As the numerical computation, there are known a method in 
which eq. (23) is directly used, and a method in which the 

15 Taylor expansion is adopted. The former is adaptive when a 
pressure change between P 0 and P is little. The later is 
suitable for an approximate expression when the pressure 
change between P 0 and P is large. When the former is 
expressed by the discrete equation where the presence is 

20 given by the present sampling time (k) and the subsequent 
sampling time (k+1), eq. (23) is expressed as follows. 
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K n { K ) n 1 < 24 > 

In eq. (24), V (k+1) denotes the estimated volume, 
and A V p (k + i) _ denotes the volume change. The Taylor 
expansion of eq. (24) is expressed as follows. 
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Here, j denotes degree of the expansion and is 
expressed by 0 • 1 • 2 * 3"-. 

(3) Result of simulation 

A result of simulation, as to the functional model 
5 of Fig. 17, which is performed with characteristic values 
of table 1, is shown in Figa. 18 and 19. Incidentally, 
with respect to the spring stiffness K s , eq. (25) is 
applied. 

10 [table 1 ] 



characteristic 


values 


of pneumatic 


cylinder 


characteristic name symbol 


unit characteristic value 


mass 


M 


[kg] 


100.0 


coefficient of 


C 


[N/(m/s)] 


1000.0 


visconsity 








esistance resistance 








coefficient of 


D 


[N/(m/s) ] 


3xl0 9 


internal damping 








ratio of specific heat 


K 




1.4 


cross-sectional area 


S 


[m 2 ] 


0.01 


cylinder volume 


V 0 


[m 3 ] 


0.001 


atmospheric pressure 


Po 


[Pa] 


0.1x10 s 



25 Fig. 18 shows time history variations of the 

respective state quantities. Fig. 19 is a view showing P-V 
and v-V diagrams representative of a relation between the 
cylinder volume V and the pressure P. In these both 
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diagrams, a solid line denotes the heat insulating state 
wherein coefficient of internal damping D = 0 , and a dotted 
line denotes the heat radiation wherein D = 3 x 10' 
[N/(m/s)]. The simulation was performed on a process in 
5 which the external step force f of 3000 [N] is applied so 

that the non- linear damping oscillation converges. In Fig. 
18, the upper side represents a response characteristic of 
an internal pressure p [pa] of a cylinder; the middle a 
response characteristic of a piston velocity v [m/ s e 

10 c ] ; and the lower side a response characteristic of a 
cylinder volume V [m 3 ] . As the common mater of those 
characteristics, it is understood that the damping 
oscillation waves are associated with a distortion owing to 
a non-linearity of a spring stiffness K [p a /m 3 ] of air. 

15 In the heat radiation state, there appears a volume change 
due to the heat shrinkage from the cylinder volume V [m 3 ] . 
' * In Fig. 19, the upper side represents a relation 

of convergent process between the cylinder volume V [m 3 ] 
and the piston velocity v [m/ sec] , in which the center 

20 of converging spiral is the stationary state. The lower 
side represents a relation between the cylinder volume V 
[m 3 ] and the internal pressure p [pa] of a cylinder, in 
which it indicates that a spring stiffness of air is of a 
non- linear characteristic. In the heat radiation state, 

2 5 it is understood that the volume change due to the heat 
shrinkage occurs regardless of pressure. 

3. Coefficient element (Oscillating mechanism) 
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As the typical example of the non-linear 
coefficient element, there is a link mechanism. 
Incorporating a mechanistic model, wherein a coefficient of 
the coefficient element (basic functional element) is given 
5 with a non-linearity, makes it possible to implement 

modeling of the non-linear coefficient. By way of example 
of the modeling, there will be discussed an oscillation 
mechanism for converting the rotary motion shown in Fig. 20 
into the oscillation motion. 

10 Fig. 20 is a view showing a structural model of an 

oscillating mechanism comprising a rotational arm A x and 
an oscillation arm A 2 . Where R denotes an arm length 
(radius); 9 a rotational angle; f an effective force on a 
coupling point; v a velocity of a coupling point; T a 

15 torque; and CO an angular velocity. The suffixes 1 and 2 
imply the rotational node and the oscillation node, 
respectively. L 0 denotes a distance between the 
respective rotary centers. While Ri is constant, R 2 is 
variable together with motion. Incidentally, both the arms 

20 are treated as rigid bodies, and inertia moment and the 
rotary resistance are neglected. 
(1) Functional model 

In modeling, a functional model, wherein a 
coefficient of the basic function element is given by the 
25 non-linear transfer factor O , is expressed by the 
following equation. 
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t 2 = *r, 

Next, there is performed modeling of a mechanistic 
model to be incorporated into the non- linear transfer 
factor <& of eq. (26). The mechanistic model of the 
5 oscillation mechanism is different in the non- linear 

transfer factor 0 according as a rotational angle for 
determining the coefficient is a rotational angle 6 x of 
the rotational arm A : or a rotational angle 6 2 of the 
oscillation arm A 2 . The non- linear transfer factors 0 a 
10 and <1>b , which are associated with the rotational angles 
By and B 2 , respectively, are in a mutually reciprocal 
relation . 

(2) Mechanistic model 1 (Oscillation model of 
rotational node input) 

15 Fig. 21 shows a functional model in which a 

mechanistic model for determining a non-linear transfer 
factor $ a with an angular velocity co 2 of an oscillation 
arm A 2 , incorporated. Incidentally, o>_ of the 
mechanistic model denotes an estimated angular velocity, 

20 and 8 2 _ denotes an estimated angle. 

Government equation and coupling condition 
equation of Fig. 21 are expressed by the following equation. 

\t 2 I U o LJ 

JL J l (2 7) 

CO = co z "I 
<I> = O fl J 



(2 6) 
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A non-linear transfer factor $ a of Fig. 21 is 



expressed by the following equation. 




5 



<D a =1- 



L 0 cos(jt - 0 2 ) 



(28) 



^ 2 -Z, o 2 sin 2 (jr-0 2 _) 



In eq. (28), the upper side denotes an estimated 
rotational angle 9 2 _ of the oscillation arm A 2 , and the 
lower side denotes the non- linear transfer factor <£> a - 



10 



Incidentally, with respect to the lower side of eq. (28), 



since all the parameters represent shapes. Fig. 21 shows a 
shape model with an octagonal frame. When the right side 
denominator of eq. (28) is 0, the oscillation arm A 2 
offers the maximum oscillation rotational angle 9 2 _max. and 
15 the non -linear transfer factor $ a is 00 . The maximum 
oscillation rotational angle 6 2 _ max is expressed by the 
following equation. 



length R 2 is expressed by the following equation, in which 
the upper side is expressed with a rotational angle 9 1 of 
the rotational arm A 1 , and the lower side is expressed 
with a rotational angle 9 2 of the oscillation arm. 




(2 9) 



20 



Incidentally, at that time, the oscillation arm 
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R 2 = ^L 0 2 +R t 2 -2L 0 R 1 cos6 1 

= L 0 cos(0 2 ) - ^R 2 -L o 2 sin 2 (0 2 ) 

(3 0) 

Equations (28) and (30) can be easily introduced 
from a geometric relation of the mechanism shown in Fig. 20, 
a load of the coupling point and a velocity vector. 

(3) Mechanistic model 2 (Oscillation model of 
oscillation arm input) 

Next, Fig. 22 shows a functional model into which 
a non- linear transfer factor ^ b , wherein an angular 
velocity CO 1 of the rotational arm A x is adopted in form 
of observation quantity, is incorporated. Incidentally, 
Fig. 22 is in a duality with Fig. 21. In Fig. 22, an 
inertia moment J and an additional function of coefficient 
C of resistance are connected to a rotational arm Ai at the 
driving side and an oscillation arm A 2 at the load side, 
respectively. From the exterior, torque T ± is applied to 
the inertia moment J to output the angular velocity C0i . 

The functional model of Fig. 22 can be expressed 
by the following domination equation and coupling condition 
equation. 
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(3 1) 

In eq. (31), first and second lines denote a state 
equation and an input and output equation, third and 
5 fourth lines denote observation equations, and fifth and 
sixth lines denote estimated observation equations for a 
rotational angle 6 X . In eq. (31), X ' B b and X BB denote 
internal storage observation quantity of an angular 
velocity co x . A mechanistic model of a non-linear transfer 
10 factor Ob of Fig. 42 is expressed by the following 
equation . 

^ R* - R,L 0 008(0, ) 

b R* +L 0 2 -2R 1 L 0 cos(d 1 ) < 32 > 

The lower side of eq. (32) denotes, similar to eq. 
(28), a shape model. 
15 (4) Result of simulation 

Next, a result of a simulation , wherein 
characteristic values of an oscillating mechanism of table 
2 are applied to the functional model of Fig. 22, is shown 
in Figs. 23 and 24. 
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[table 2) 



characteristic values of oscillating mechanism 



characteristic name 


symbol 


unit characteristic value 


inertia moment 


M 


[N/m 2 ] 


0.1 


coefficient of 


C 


[N/(m/s) ] 


0.3 


resistance 








radial of rotational 


Ri 


[m] 


0.02 


arm 








distance between 


L 0 


[m] 


0.05 


centers of rotation 









Fig. 23 shows a response characteristic where step 
torque Tj of 0. 0 5 [Nm] is applied. In Fig. 23, the 

15 upper side denotes an angular velocity (D x [r ad/s e c] of 
the rotational arm, the middle denotes an angular velocity 
w 2 [r ad/s e c] of the oscillating node, and the lower 
side denotes a torque T 2 [Nm] of the oscillating node. 
The respective state quantities are associated with a 

20 waveform having a distortion by a non- linear coefficient 
$ b . 

Fig. 24 is a diagram of a transmission 
characteristic to an angle of a non-linear coefficient <& b . 
A distortion of the state quantity of Fig. 23 is determined 
25 by the transmission characteristic shown by this diagram. 
The arrow in the figure denotes a change direction of the 
non- linear transfer factor ^ h when the rotational arm Ai 
rotates counterclockwise. 

4. Example of modeling of friction 
30 Friction is a typical example of a non-linear side 
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load. Fig. 25 shows a mechanistic model of friction in 
which a frictional force F v is effected on mass M. Fig. 
25 shows a structural model of mass M which moves on a 
plane at velocity Vj* v 2 while external force fx* f 2 is 
5 effected by a frictional force F v . The frictional force 
F v consists of a stationary frictional force F s which is 
to be a constant and a dynamic frictional force F M which 
is to be a function of velocity. 

(1) Functional model 

10 Fig. 26 shows a functional model in which a side 

load F v generating a frictional force is applied to mass M 
of the potential storage element. A mechanistic model 
generating a frictional force is incorporated into the side 
load F v . The mechanistic model shown in Fig. 26 

15 generates a frictional force F f generated from velocity 
(internal state quantity) of mass M and an estimated 
observation quantity of an acceleration force f a • In the 
output velocity v 2 side, coefficient C of resistance 
generating reaction f 2 is added in form of a load. 

20 A government equation and a coupling condition 

equation for the functional model, which are determined 
from Fig. 26, are expressed by the following equation. 
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(3 3) 



In eq. (33), from the top in order, the first line 
denotes a government equation; the second line an input and 
output equation; and third and fourth lines estimated 
observation equations for reaction f 2 of coefficient of 
resistance and an acceleration force f a_ . In the 
coupling condition equation of the lower side of eq. (33), 
the first line denotes quantity of an observation of 
velocity; the second an estimated acceleration force to be 
applied to mass; and the third a frictional force. 

In the mechanistic model of Fig. 26, the lower 
side shows a model for generating a static frictional force 
F s ; the upper center a model for generating a kinetic 
frictional force F M ; the left end a model for selectively 
output ting quantity F f of substitution between the static 
frictional force F s and the kinetic frictional force F M . 
At the right end, a characteristic of the kinetic 
frictional force F M depending on velocity is incorporated 
on a graphic representation basis . In the event that the 
kinetic frictional force F M is a constant value 
regardless of the velocity v, no graph is needed. The 
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graph of the frictlonal characteristic shows that the 
frictional force varies from the high static frictional 
force to the low kinetic frictional force, and it is 
possible to reproduce stitching phenomenon on a frictional 
5 surface and an oscillation phenomenon wherein the 

frictional change of the time -lapse area is regarded as the 
negative resistance. 

(2) Mechanistic model (frictional model) 
Selection between the static frictional force F s 

10 and the kinetic frictional force F M is performed in 
accordance with a determination condition S w of the 
velocity v. In selection, on the stop state of S w =0, 
the static frictional force F s is selected, and on the 
moving state of S w = 1 , the kinetic frictional force F M is 

15 selected. The mechanistic model is expressed by the 
following equation. 

if (v s =0) then (5 W = 0) else (5 W =1) 

= ^Mi^WVJ. + ^Si^WV_Q 

On the second line of eq. (34), F Mi denotes a 
20 kinetic frictional force; F M and F sl a static frictional 

force; and F s and v s velocity. Suffixes 0 and 1 of switch 
variables denote an NO switch side and an NC switch side of 
a selection switch, respectively. The determination 
condition equation at the upper side of eq. (34) is 
2 5 expressed by the following equation in accordance with a 
determination method of converging it to zero. 



(3 4) 
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if (|Vs(jfc)|<|v 5W -v 5(Jt _i)|) f/iew OW =0) else (5 W = 

(3 5) 

In eq. (35), v s(k) denotes the present velocity; v 
5 s{k-i) "the velocity before one sampling period; and Av a 
velocity change width determined from these. A 
determination result of condition determination S w is the 
same as eq. ( 34 ) . 

Next , there will be introduced a mathematical 

10 model of generating the static frictional force F s . The 

static frictional force F s is associated with two behaviors. 
One of the two behaviors is concerned with such a matter 
that in the event that the acceleration force f s is not 
more than the static frictional force F s , the static 

15 frictional force F s balanced with the static frictional 

force F s is applied to mass M in form of reaction force so 
that the stop state is kept. Another is concerned with 
such a matter that when the acceleration force f s , which 
is not less than the static frictional force F s , is 

20 applied, the kinetic frictional force F M is immediately 

selected so that mass M is activated. At that time, since 
the mass is in the stop state, a direction (code) of the 
kinetic frictional force F M is the same direction (code) as 
the acceleration force f s , and the side load F v offers a 

25 negative value to be a reaction force with respect to the 
mass M . Selection between both the behaviors is performed 
in accordance with an operation of a condition decision S WY 



96 



for difference Y f = I f s I — F s0 between the acceleration 
force f s and the side load F s0 . In selection, the balance 
state of the stop is expressed by S WY =0 so that the 
acceleration force f s side is selected, and the moving 
instant is expressed by S WY = 1 so that the side load F so 
side is selected. This mechanistic model is expressed by 
the following equation. 



second line of eq. (36) is divided by I f s I is that a 
direction (code) of friction is determined. Finally, 
modeling is performed on a mechanistic model for generating 
the kinetic friction F M . The kinetic frictional force F M 
is generated in such a manner that a code of velocity v v 
is added to the side load F M0 of the kinetic frictional 
force . 

This mechanistic model is expressed by the 
following equation. 



if i(|/s|-^o)>0} then (S w =l) else (5 W =0) 




(3 6) 



The reason why the acceleration force f s on the 



F, 



M 




(3 7) 



The static frictional force F so of the eg. (36) 
and the kinetic frictional force F M0 of eq. (37) need a 
software to be determined through a graphic retrieval by 
velocity v or a mathematical model . 

The above-mentioned mechanistic models are coupled 
by the . following coupling condition equation. 



1 SO — MO " 

f=fs 



(3 8) 



(3) Result of simulation 
r - Figure 27 is a view showing a result of a 

simulation which is performed on the functional model with 
the characteristic values shown in table 3. Incidentally, 
a decision of velocity v = 0 was performed in accordance 
with eq. (35), and the kinetic frictional force F M0 was 
constant . 



20 
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[table 3] 

characteristic values of friction 



characteristic name symbol unit characteristic value 

5 mass M [kg ] 0.1 

coefficient of C [N/(m/s>] 0.5 

resistance 

static , frictional force F S o EN] 0.5 

kinetic frictional force F M0 [nj 0 . 1 

10 

In the simulation, the external force f rising at 
0.5 [N/sec] per unit time is applied up to 1 [N] , and the 
external force up to -1 [N] is applied at the same 
variation rate from 0.4 [sec] and the external force is 

15 removed at 10 [sec]. In Fig. 27, the first denotes an 
external force fi [N] ; the second a velocity v L ; the 
third a frictional force F v [N] applied in form of a side 
load in the third; and the fourth a reaction f 2 of 
resistance coefficient C . Incidentally, the frictional 

20 force F v and the reaction f 2 denote quantity of 
observation . 

In Fig. 27, applying the external force fi slowly 
causes the static frictional force F v to occur to cancel 
the external force f 1 so that mass M does not move. From 
25 the point of 1 [s e c] , the external force f 1 exceeds the 
frictional force F v so that the mass M starts to move. 
Thus, a certain kinetic frictional force F v generates and 
it reaches the stationary state. From the point of 4 [s e 
c] , the external force fi slowly downs, and from 6 [s e 
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c] or so, again the static frictional force F v generates 
so that the mass M temporarily. Further, when the 
external force f x is applied in the opposite direction, 
the mass M starts to move in the opposite direction so that 
5 the kinetic frictional force F v generates. Finally, the 
external force f x is removed, the mass M is lowered in 
velocity V ± to be stopped. In stopping, velocity due to 
the force of inertia of the mass M appears in form of 
reaction f 2 via the coefficient C of resistance, while the 
10 static frictional force F v generates, and the mass M stops 
after a while. After the stop, the frictional force F v 
disappears . 

5. Potential storage element (variable moment of 

inertia) 

15 As a mechanism in which moment of inertia offers a 

non-linearity, there is a mechanical element such as a 
centrifugal governor and a centrifugal clutch. Fig. 28 is 
a view showing such a mechanical element by way of example . 

In Fig. 28, torque T x , where a direction shown in 

20 the figure is defined as the forward rotation, is effected 
to generate an angular velocity co x . The rotational system 
shown in Fig. 28 is a structural model into which a system 
where a movable mass M m receives a centrifugal force to be 
translated is incorporated. The movable mass M m rotates on 

25 the circumference of the radius R m , so that the 

centrifugal force F e acts on the movable mass M m to move 
outside at velocity V e . According to this rotational 
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system, when the system rotates, upon receipt of torque T\ , 
at the angular velocity a>i , velocity of the movable mass 
M m in the circumference direction is given by V m , and 
force of inertia in the circumference direction is given by 
5 F m . The lower limit (stop position) of the rotational 

radius R m of the movable mass M m is expressed by R 0 • In 
the translational system, a rotation of the movable mass 
M m causes the centrifugal force F e to be generated, so that 
the movable mass M m moves by L ra at velocity V e to offer a 

10 rotational radius R m , while receiving a resistance force 
of a coefficient C m of viscosity resistance against 
stiffness K m of a spring. Variation of this rotational 
radius R m causes moment J of inertia of a rotational 
system to be varied. When the movable mass M ra moves to 

15 the radius R s , it contacts with a stopper of stiffness K s 
and stops. Where J 0 denotes moment of inertia of a 
suspension frame excepting the movable mass M m , stiffness 
K m of the spring has no mass. Incidentally, stiffness of 
stop position R 0 is not considered. 

20 (1) Functional model 

An aspect of the mechanism of a variable moment of 
inertia resides in the point that mass M causes the moment 
J of inertia to be changed by the centrifugal force F e . 
Consequently, the basic structure of the functional model 

25 has a mechanistic model wherein moment J m of inertia to 

which variable correspondence of a moving distance L m of 
the movable mass M m moved by the centrifugal force F e due 
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to the observation quantity of the estimated angular 
velocity CO !_ of the rotational system is added, is 
substituted for moment J of inertia of a potential storage 
element (a basic functional element). And this mechanistic 
model includes a functional model of the movable mass M m 
moving by the centrifugal force F e . Fig. 29 shows 
functional and mechanistic models for a rotational shaft 
representative of this relation. In Fig. 29, a coefficient 
C of viscosity resistance is added to the basic element of 
the potential storage element. Further, in Fig. 29, an 
output potential system of the mass M m includes a potential 
quantity of moving velocity V m and a storage potential 
quantity of rotational radius R m . 

Government equation and coupling condition of the 
functional model shown in Fig. 29 are expressed by the 
following equation. 



0 " 




-j -c r 


CO, 




0-10 



CO j = co 1 

J =J m 

(3 9) 

(2) Mechanistic model (centrifugal force model) 
The mechanistic model comprises a rotation to 
translation transformation for coupling a rotational system 
with a translational system, a translational motion of a 
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movable mass Mm, and a function of a stopper of the movable 
mass M m , which determines the rotational radius R m of the 
movable mass M m . 

First, the mechanistic model of the rotation to 
5 translation transformation is expressed by the following 
mathematical model. 

J m =M m R m 2 + J 0 \ 

In eq. (40), the upper side denotes centrifugal 
10 force F e of the movable mass M ra by the rotational system, 
and the lower side denotes moment J m of inertia varying in 
accordance with the centrifugal force. 

The functional model, wherein the movable mass M m 
is subjected to the translational motion in accordance with 
15 the centrifugal force F e , is the following government 
equation. 
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In eq. (41), the upper side denotes the government 
equation of the translational motion in which K ra and K s in 
Fig. 29 are unified through an addition coupling of the 
integration quantity. And the lower side denotes an 
5 estimated observation quantity in which the storage 

potential quantity of the movable mass M m is expressed by 
the moving distance L ra . Incidentally, t smp is a sampling, 
period. 

From eq. (41), the estimated observation quantity 
10 of the rotational radius' R m is expressed by the following 
equation . 

R m = ^0 + L m ( 42) 

Finally, the stopper for limiting movement of the 
movable mass M m is expressed by the following equation. 

15 

L s =R s ~ R m 

if {L s i> 0) then (S m = 1) else (S m = 0) 

In eq. (43), condition determination at the time 
when it is in contact with the stopper is given by the non- 
20 contact state Sws^O, and the contact state S ws = 1 . 

The coupling condition of assembling these 
mechanistic models in form of the whole mechanistic model 
is expressed by the following equation. 



(4 3) 
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v e =v m =v s 

(3) Result of simulation 

A result of the simulation performed on the 
5 functional model of Fig. 29 using the characteristic values 
shown in table 4 is shown in Fig. 30. 

[table 4] 

characteristic values of variable moment of inertia 



ffl 





10 


characteristic name 


symbol 


unit characteristic value 






rotational system 












inertia moment 


Jo 


[N/sec 2 /rad] 


0.01 






(suspension frame) 








o 




coefficient of 


C 


[N/(m/sec) ] 


2.5 




15 


viscosity resistance 












translational system 












movable mass 


M m 


[N/sec 2 /m] 


0.2 






coefficient of 


Cm 


[N/(m/sec) ] 


10.0 




20 


viscosity resistance 












stiffness 


K m 


[N/m] 


2.0 






stopper stiffness 


K s 


[N/m] 


1.0 x 10 4 






rotational radius 


Ro 


[m] 


0.1 






lower limit 










25 


rotational radius 


R s 


[m] 


0.5 



upper limit 



In Fig. 30, according to the simulation, first, 
30 input torque Tj.= 2 0 [N] is applied, after 0. 4 [sec] 
torque l\=-2 0 [N] is applied in the opposite direction 
to rotate it reversely, and after 0. 9 [s e c] , torque T x 



(44) 
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= 2 0 [N] is applied in forward direction, and about 1 . 4 5 

[sec] it contacts with the stopper. In Fig. 30, 
regarding results, from the top in turn the first denotes 
an angular velocity [r a d/s e c] of the rotational 
5 system; the second a centrifugal force F e [N] acting on 

the movable mass M m ; the third a moving velocity V e [m/ S 
ec] ; the fourth spring and stopper stiffness force F s 
[N] , and the fifth a rotational radius R m [m] of the 
movable mass M m . 

10 In the simulation result, selection of input 

torque T x causes the centrifugal force F e of the movable 
mass M ra to be 0 in an instant and again return. While the 
moving velocity V e is lowered in an instant of selection, 
it returns so as to continue to move. This is independent 

15 of the rotational direction. When the rotational radius 
R m offers 0 . 5 [m] , the movable mass M m contacts with 
tire stopper. In the contact state, repulsion due to a 
collision causes damping oscillation to be generated on 
moving velocity V e and force F s of K s of a stopper 

20 stiffness. At that time the damping oscillation wave is 
flat on the lower side. It is understood that this is a 
repulsion force in an instant that the movable mass M m 
contacts with the stopper. 

Figure 31 is a principle explanatory view of a 

25 second non-linear characteristic reproducing apparatus 
according to an embodiment of the present invention. 
Figure 31 is a view for easily understanding the principle 
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of the second non-linear characteristic reproducing 
apparatus according to the present invention, and thus the 
second non-linear characteristic reproducing apparatus 
according to the present invention is not confined to Fig. 
5 31 and the associated explanation. 

Figure 31 shows a model structure in the event 
that the concept of the first non- linear characteristic 
reproducing apparatus according to the present invention is 
applied to the whole system. According to the conventional 

10 modeling of a non-linearity, an output state quantity such 
as characteristics and coefficients is controlled by 
functional values of input state quantity such as an 
exponent, a logarithm and a square. Consequently, it is 
difficult to solve a plurality of non -linear physical" 

15 phenomena through a mathematical model integration. 

However, the use of this method makes it possible to 
execute computation regarding characteristic values of the 
state quantity transformation sections Ri, R 2 and R 3 as 
being linear, at the moment of execution of the whole 

20 system model, since the parameters (characteristics, 

coefficients and attached load) for producing non-linear 
physical phenomena are set up through determination of 
prediction values of the state quantity transformation 
sections R lt R 2 and R 3 prior to execution of the 

25 computation. 

(Explanation of examples) 
6 . Geneva mechanism 
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(1) Basic function 

As one of the mechanical elements to which an 
oscillating link mechanism is applied, there is a Geneva 
mechanism for performing a sampling action for a feed 
5 controlling device or the like. The Geneva mechanism 

comprises a master section and a slave section. The master 
section is associated with a rotational arm of the 
oscillating link mechanism. The slave section is provided 
with a groove opening on the circumference side in 

10 association with an oscillating arm. A connecting shaft of 
the master section is engaged with the groove of the slave 
section so that rotation of the master section causes the 
slave section to intermittently perform a stop and a 
rotational motion. Fig. 32 shows the Geneva mechanism. 

15 The Geneva mechanism of Fig. 32 is given with 

torque T p • T s and angular velocity co p • <x> s on the master 
section and the slave section, respectively, and is 
provided with moment of inertia J p * J s and coefficient of 
viscosity resistance C p • C s on the master section and the 

20 slave section, respectively. In Fig. 32, L z denotes a 

distance between the nodes, and R p denotes a rotational 
radius of a connecting shaft. In application of the 
rotational torque T p to the master section, when a 
rotational angle j3 w of the connecting shaft of the master 

25 section is placed at the position of a connecting angle — 
a of the connecting groove of the slave section, both the 
elements are engaged with one another to drive the slave 
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section. And when the slave section rotates to the 
connecting angle 2 a , the slave section is disconnected and 
stopped. In the stopped state, a semicircular shape of 
cut-out of the slave section and a circular member of the 
5 master section are engaged with one another, so that the 

slave section cannot be rotated. Incidentally, the number 
of intermittent operations of the slave section to one 
rotation of the master section is determined by the number 
N of connecting grooves provided on the slave section. 

10 (2) Functional model (Geneva model) 

Fig. 33 shows function and mechanistic models in 
which the structure of Fig. 32 is subjected to modeling. 
In Fig. 33, T p and co p denote drive torque and rotational 
angular velocity of the master section, respectively, T s 

15 and co B denote drive torque and rotational angular velocity 
of the slave section, respectively, and an observation 
quantity 6 s denotes a rotational angle of the slave 
section. In the functional model of Fig. 33, the master 
section is provided with moment of inertia J p and 

20 coefficient C p of viscosity resistance, and the slave 
section is provided with moment of inertia J s and 
coefficient C s of viscosity resistance. The angular 
velocity C0 P of the master section is applied via a transfer 
factor 0 to the non-linear stiffness K, and then to the 

25 slave section through conversion to torque. The torque is 
returned via the transfer factor <f> to the master section 
in form of a load torque. The government equation of the 
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functional model of Fig. 33 is expressed by the following 
equation. In Fig. 33, x denotes an internal state quantity. 
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(45) 

The mechanistic model of Fig. 33 is of a 
structural non-linearity. Consequently, the rotational 
angle |3 w of the master section is determined in form of an 
estimated observation quantity, and included in an output 
potential quantity of moment of inertia J p . The rotational 
angle 3 W is expressed by the following equation. 



dt 



(4 6) 



The coupling condition equation for causing the 
mechanistic model to be nesting in the functional model is 
expressed by the following equation. 
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K 



(4 7) 



5, 



(3) Mechanistic model (intermittent motion model) 



The mechanistic model incorporated into the model 



5 



of Fig. 33 includes a function of a coupling condition S 



of determining a rotational angle B w of the master section 
to decide a connection with the slave section, a function 
of determining a connecting stiffness K z in the connection, 
and a function of determining transfer factor <t> z between 
10 the master section and the slave section. These 

subordinate mechanisms are subjected to a horizontal 
expansion in accordance with the following coupling 
condition equation. 



15 




(4 8) 



(a) Mechanistic model (Connection deciding model) 



From the coupling condition equation (47), the 



estimated rotational angle j3 w _ of eq. (46) is inputted to 



20 



the mechanistic model in form of the estimated rotational 



angle 0 B _. The estimated rotational angle 0 B _ is able to 
be transformed in an angular range of 0 ~ 2 % in form of 
the rotational angle of the slave section. The j3 is 



inputted into the respective mechanistic models in 
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accordance with eq. 



(48) . 



= 2jz 



2jz 



-mt 



2ji 



L + 0.5 



(4 9) 



In eq. (49), i n t (a) denotes a function of 
cutting off a decimal part of value 'a' and deriving an 
integer part. In the structure model of Fig. 32, a 
connection angle a for connection of the slave section to 
rotation of the master section is determined in accordance 
with the number N of connecting grooves of the slave 
section to the connecting shaft of the master section. The 
equation is as follows . 
Jt 

]v~ < 50 > 

In eq. (50), The connection angle a p of the master 
section to the connection angle a of the slave section can 
be expressed by the following equation from the geometric 
structure of Fig. 32. 



a = 



-i \Lz 



( 



sin 2 (a) + 



-cos(a) 



-^sin(2a) 



(5 1) 

In eq. (51), in the event that the rotational 
angle j3 is within a range ± a p of the connection angle. 



112 



the master section and the slave section are connected with 
each other. Consequently, the condition decision S ws of 
determining the connection state is expressed by the 
following equation. 

if (P>a p ) then ( S ws = 0) else (S ws = 1) 

( 5 2 ) 

In eq. (52), S ws = 1 denotes a connecting state; S 
ws = 0 a non-connecting state. 

(b) Mechanistic model (transfer factor model) 

Next, the transfer factor <f> z between the master 
section and the slave section in the connecting state is 
associated with the equation of the link mechanism 
identical with the oscillating mechanism, and thus is 
determined by a shaft -to -shaft distance L z of the master 
section and the slave section, a rotational radius R p of 
the connecting shaft of the master section, and the 
rotational angle j3 of the master section. The 
mathematical model is expressed by the following equation. 



R -R„L z cos(/3j 

P P * v/ <t>/ 

q> - - -- '- ~ ' " (5 3) 



R; +L/ -2R p L z cos(^) 



In eq. (53), the transfer factor cp z shows that at 
the time of non-connection the switch variable S w shown in 
Fig. 33 becomes 0 so that the master section and the slave 
section are disconnected from one another. 
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{ c ) Mechanistic model (transfer stiffness model) 
The non-linear rigidity K between the master 
section and the slave section is a combination of stiffness 
of the master section, stiffness of the slave section and 
5 stiffness of the contact. The non-linear stiffness K in 

connection offers stiffness via the connecting shaft of the 
master section and the connecting groove of the slave 
section, and can be represented by a function of the 
rotational angle j3 w of the master section. In the non- 
10 connection, the non- linear stiffness K offers stiffness of 
the top of the connecting groove of the slave section, the 
circular configuration of the master section and the 
contact point. When stiffness is large, it is possible to 
deal with the non-linear stiffness K as constant. 
15 (4) Result of simulation 

Figures 34 and 35 show results of a simulation in 
which the simulation for functional and mechanistic models 
is performed using characteristics of table 5. 



25 
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In Fig. 34, the result is represented as a time 
history. In Fig. 35, the horizontal axis represents the 
rotational angle (3 w of the master section. In Fig. 34, 
from the top in turn, co p and co s [r e d/ sec] denote 
5 angular velocities of the master section and the slave 

section, respectively; x K [Nm] a torque of the non- linear 
stiffness K; and S w a condition determination. In Fig. 35, 
the horizontal axis denotes the rotational angle 0 W [r a 
d] of the master section. With respect to the vertical 

10 axis, at the upper side, co p and co s [r e d/ sec] denote 

the angular velocity in a similar fashion to Fig. 34, and 
at the lower side, 9 S [r ad] denotes the rotational angle. 
The simulation is associated with a result in the event 
that 2. 7 [Nm] is applied to the drive torque T p of the 

15 master section on a step basis, and the load torque T s is 
not applied to the slave section. 

In Fig. 34, when the connection is established 
through the switch variable S w = 1 , the master section 
drives the slave section, so that the angular velocity co p 

20 of the master section is remarkably lowered because the 

slave section becomes the load. At that time, the angular 
velocity of the slave section rapidly rises. However, 
lowering of the master section in the angular velocity 
involves lowering of the angular velocity C0 S on the slave 

25 section. And the torque x K of the non-linear stiffness K 
generates slightly an oscillation owing to an impact of the 
connection. While the slave section stops on the condition 
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of condition decision S w =0 at the time of the connection 
termination, a slight oscillation is generated on the 
torque x K owing to an impact of the stop. 

From Fig. 35, it is understood that the rotational 
5 angular 6 s of the slave section intermittently operates in 
unit of k/2 every a rotation ( 2 7t) of the master section, 
and in the non-connection the rotational angle is 
maintained. 

7. Modeling of a residue warning lamp 

10 According to a complicated system of general 

products and parts, it is permitted that a slow transient 
response is regarded as a non- linear in which a stationary 
characteristic (static characteristic) is slowly changed, 
and is replaced approximately by a mechanistic model to be 

15 incorporated into a linear functional model of a high speed 
response. As an example, there is raised a residue warning 
lamp utilizing a contactless switch of a negative 
characteristic resistance. 

(1) Principle of detection of liquid level 

20 Fig. 36 shows a part constitution of a residue 

warning lamp and its connection. The operation of the 
residue warning lamp is associated with a very simple 
system in which when a height of a liquid level in a tank 
shown in Fig. 36 is above L 0 , the warning lamp turns off, 

25 and when it is below L Q , the warning lamp turns on. The 
used thermistor is a device which varies in a resistance 
value in accordance with temperature. Such a device is a 
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general non- linear resistance element which is widely used, 
in an electronic circuit and a temperature measurement 
system. This system is required in a functional aspect of 
preventing a flicker of the warning lamp due to shaking of 
5 a liquid level caused by start and stop of a car and 
jolting of a car, and an aspect of quality of high 
reliability along with a low cost. These requirements are 
satisfied by utilization of the non-linear characteristic 
of the thermistor. 

10 Figure 37 is a circuit diagram of an electric 

circuit of a system of the residue warning lamp. 

In Fig. 37, R B denotes an internal resistance of a 
battery; E 0 an internal electromotive force of the battery; 
R L a resistance of the warning lamp; and R T a resistance of 

15 the thermistor. These elements are connected in series. 
With respect to state quantities, V B denotes a power 
source voltage; V L a terminal voltage of the warning lamp; 
V T a terminal voltage of the thermistor; and I B a current 
common to these three parts . 

20 Prior to modeling of the residue warning lamp, 

there will be simply described an operation of a 
contactless switch in which a thermistor resistance R T 
operates a warning lamp resistance R L of the load. Figure 
38 shows voltage -current characteristics of R T and R L and 

25 the relation. 

In Fig. 38, it is considered that voltage V L of 
resistance R L is a voltage drop from the power supply 
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voltage V B , and there is shown voltage -current 
characteristics in the opposite direction in which the 
original point is shifted to point V B on the voltage axis. 
A cross point of R T and R L offers current I B conducting in 
5 the circuit. The left side of the cross point represents a 
voltage V T to be applied to the thermistor, and the right 
side to V B represents a voltage V L to be applied to the 
warning lamp. The thermistor resistance R T has a 
resistance characteristic that it is larger at a low 

10 temperature and smaller at a high temperature. Thus, the 
thermistor serves as a switch wherein the cross point of 
Rt_off . which is larger in resistance value, offers an off- 
state, and the cross point of R t _on , which is smaller in 
resistance value, offers an on-state. Both the resistors 

15 are connected in series with one another. Thus, in the 

off-state, the current I b _off denotes a leakage current; V 
t_off a cut-off voltage by the thermistor; and V L O ff a turn- 
off voltage of the warning lamp. In a similar fashion, in 
the on-state, a current I B _ ON denotes a turn-on current; V 

20 t _on a voltage drop of the thermistor; and V L O n a turn-on 
voltage of the warning lamp. 

In the thermistor, as shown in Fig. 37, a current 
I B always conducts through the warning lamp, so that 
consumed power W T of V T I B or I b 2 Rt causes a self -heat 

25 generation to occur. This involves a temperature rise. In 
the event that the liquid level is high, the thermistor is 
in liquid, and thus the generated heat is cooled by the 
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surrounding liquid, so that the thermistor resistance R T 
becomes large and thereby offering the off -state in which a 
circuit current I B is small . When the liquid level is 
lowered, so that the thermistor appears, the thermistor is 
5 not cooled. Accordingly, the temperature rises, so that 
the thermistor resistance R T is lowered. Thus, an 
increment of the circuit current I B causes the consumed 
power W T to be increased to generate a heat. Further, 
lowering of R T causes an increment of I B whereby a 

10 thermal runaway occurs. In the thermal runaway, owing to 
the suppression of the circuit current I B by the 
resistance R L of the warning lamp, the terminal voltage 
V T of the thermistor becomes smaller according to 
transition to the on-state, so that the consumed power W T 

15 is lowered and whereby the on- state is maintained at the 

balancing point of self -heat generation and heat radiation 
of the thermistor. Here, the resistance R L of the warning 
lamp has a non-linearity of the positive characteristic 
that when the voltage is low, the resistance is small, and 

20 when the voltage is high, the resistance is large, contrary 
to the thermistor as shown in Fig. 38. In the residue 
warning lamp, the non-linearity of both the static 
characteristic and the negative characteristic is utilized 
to implement the off -state in which a leakage current is 

25 less and the on-state in which a voltage drop is less. 

With respect to the flicker of the display as mentioned 
above, the flicker is prevented by utilizing the response 
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delay to the temperature of the thermistor transiting 
between ON/OF F . 

(2) Functional model 

Fig. 39 shows a functional model in which the 
above-mentioned contents are rearranged to implement a 
device modeling and a mechanistic model is incorporated. 

In Fig. 39, a battery is a current supply for 
supplying a current through an internal resistance R B in 
accordance with a difference between a source voltage E 0 
and a terminal voltage V B . A resistance R L of the 
warning lamp has a mechanism depending on a voltage, and a 
thermistor resistance R T has a mechanism depending on a 
consumed power. The temperature rise T x due to the self- 
heat generation is included in the output potential 
quantity of the thermistor. Details of the respective 
mechanistic models will be described later. 
" ' A mathematical model is derived from Fig. 39 as 

follows . 

The mathematical model of the battery is expressed 
by the following equation. 



i b -4-( e o- v *) 



r> v - (54) 



The respective terminal voltages of a series 
2 5 circuit for the warning lamp and the thermistor is 
expressed by the following equation. 
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V T = R T I B 
V L = RJ B 



(55) 



(5 6) 



A result of substitution of eqs . (54), (55) and 
(56) for one another and rearrangement is as follows. 
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(57) 



(5 8) 



From eqs. (57) and (58), the government equation 
of the residue warning lamp is expressed by the following 
equation . 



V 




' g 0 e; 















[1] 



(59) 



where G 0 denotes a overall conductance of the 
whole system expressed by the following equation. 



1 

R b+ R T+ R L 



(6 0) 
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Equation (59) is a stationary characteristic 
constituted of an observation equation having no state 
equation and input /output equation. 

(3) Mechanistic model 

(a) Mechanistic model of a warning lamp 
As the warning lamp, there is used an automotive 
incandescent lamp where rated voltage is 12 [V] and rated 
power is 1.4 [W] . Table 6 shows specifications. 

[table 6] 

specification of warning lamp 

characteristic name symbol unit characteristic value 
warning lamp : 

power W L [W] 1.4 

voltage V L [V] 14.0 

current I B [A] 0.1 ± 0.01 

flux of light $l [lm] 90.4 

size : 

length L [mm] 19.5 
width W [mm] 11.3 
height H [mm] 11.0 

It is known that voltage-current characteristic is 
expressed by the following equation. 



(61) 
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In eq. (61), V s denotes a reference voltage to set 
up characteristics, I s denotes a reference current in the 
reference voltage V s . From eq. (61), a non-linear 
resistance of warning and the like can be introduced as 
5 follows . 



i. vav* (62) 



A non-linear resistance, in which eq. (62) is 
10 rearranged in accordance with Ohm's law, is expressed by 
the following equation. 



Rj = R.n I TZ~ I (63) 



^, 

Where R s denotes a reference resistance in the 
15 reference voltage V s , and is expressed by the following 
equation . 

v s 

R s = ~T~ (64) 
1 s 

Figure 40 shows a mechanical resistance of a non- 
20 linear resistance R L of eq. (64). 

Now let us determine the resistance R L of the 
warning lamp and the multiplier factor n from eq. (63) and 
the measured values of the voltage and the current, 
respectively. Figure 41 shows a result of measurement 
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where a dot line denotes measured values , and a solid line 
denotes a computed value. From this result, the multiplier 
factor n determined by least squares method is n^O, 5 5, 
and the reference current, where the reference voltage V s = 
1 4 [V] , is I s =0. 1 0 2 [A] . 

(b) Mechanistic model of negative characteristic 
resistance element 

Fig. 42 is a perspective view of a thermistor. 
Fig. 43 is a view showing a mechanistic model of the 
thermistor. 

The thermistor shown in Fig. 42 is housed in a 
cylindrical metallic case having three apertures through 
which liquid flows. 

In Fig. 43, a necessary mathematical model is" 
introduced from the characteristic equation to implement 
modeling. The non-linear resistance characteristic of the 
thermistor is expressed by the following equation depending 
on a temperature . 

R T = R 0 qx^b t ^~y^ ( 65) 

In eq. (65), T 0 denotes a reference temperature; B 
T a B constant to determine a thermistor characteristic; 
R 0 a reference resistance in the reference temperature T 
o ; and T x a used temperature. 

The rising temperature T a due to the self-heat 
generation of the thermistor is expressed by the following 
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equation with heat dissipation factor 6 T and consumed power 
W T . 

T u ~ + T o (66) 

From the functional model of Fig. 39, the consumed 
power W T of the thermistor is expressed by the following 
equation. 

W T = I B V T 

r 2 p (6 7) 

While it is possible to determine a resistance 
value R T from the temperature rise due to the consumed 
power W T of the thermistor in accordance with eqs . (6 5) , 

(66) and (67) , these equations are of a stationary 
state after the temperature rise. For this reason, from 
the heat time constant (a dry state) of the thermistor, 
established is modeling of the response characteristic 
associated with the situation that the thermistor moves in 
liquid near the practical use state and in the air. This 
state is concerned with two conditions one of which is a 
state wherein the thermistor is dried from the wet, and 
another is the opposite state. In view of this situation, 
a transient model of the secondary delay is applied 
considering delay of rising up. As an equation of the 
secondary delay, the following discrete equation, which is 
normalized in range of 0 - 1 in a similar fashion to that 
of the temperature rise model as aforementioned, is used. 
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U(k) 



(6 8) 



Equation (68) is a secondary response 



characteristic equation in which two primary delay of 
discrete equations interfere with one another. In eq. (68), 
P a - P b denote discrete system inherent values having values 
of the range of 0 ~ 1 , and G ab denotes a factor of 
distributing two primary delay of inputs. Equation is as 
follows . 



In eq. (69), t smp denotes a sampling period, and T 
a and r b denote a time constant. 

Between a case where the thermistor is picked up 
from liquid to the air and a case where the thermistor is 
entered from the air into liquid, the heat radiation is 
different owing to the difference in the way of wetness, 
and thus, of course, also the response characteristic is 



P a 



exp{-^ /tJ 



*a + *b 



(6 9) 
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different. In view of the foregoing, <5 T of eq. (66) and X 
a , r b of eq. (69) are changed over between the liquid and 
the air. A switch element for selecting between the two 
states is denoted by S w . <5 T , r a , and r b , wherein the 
liquid and the air are identified by w a t and dry of a 
subscript suffix, respectively, are expressed by the 
following equation. 



^ a ~ ^a_wat^W_0 ^ a _dry^W _1 
T b = T b_wat^W_0 +T b_dry^W_l 



T_wat^W_0 _dry^W _1 



d T =S 



(7 0) 



In the switch element of eq. (70), the liquid is 
denoted by S w _ 0 =land S w _!= 0 , and the air is denoted by S 
wjq— 0 and S w _i = 1 . 

(4) Result of simulation 

A simulation of a residue warning lamp is 
performed in such a manner that the characteristic values 
of table 7 are applied to equation of the functional and 
mechanistic model. Figures 44 and 45 show the result of 
such a simulation. 
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Figure 44 is a view showing a step response in the 
event that the thermistor is put into the liquid and is 
picked up therefrom into the air. In Fig. 44, from the top 
in order, S w denotes a condition determination; V T • V L 
5 terminal voltages of the thermistor and warning lamp; I B a 
circuit current; and T x • W T temperature and consumed power 
of the thermistor. 

Figure 45 shows a response result of a 
reproduction wherein the thermistor is covered by liquid at 

10 a period of 2 [sec] owing to a vibration of the liquid 
level, which is expressed on a similar basis to Fig. 44. 
In Fig. 45, a time ratio of the dry in one period is varied 
at 0 [%] ~5 0 [%] ~1 0 0 [%] and, it is continued by 
3 0 [ s e c ] in case of 5 0 [%] , . 

15 Figure 46 is a principle explanatory view of a 

third non-linear characteristic reproducing apparatus 
according to an embodiment of the present invention. 
Figure 4 6 is a view for easily understanding the principle 
of the third non-linear characteristic reproducing 

20 apparatus according to the present invention, and thus the 
third non-linear characteristic reproducing apparatus 
according to the present invention is not confined to Fig. 
46 and the associated explanation. 

The apparatus of Fig. 46 comprises a state 

25 quantity selecting section (switch element) S w for 

performing connection and disconnection of a state quantity 
between input state quantity and output state quantity of 
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the same physical quantity, and a logical decision section 
for generating a signal to disconnect the state quantity 
selecting section. The state quantity selecting section is 
given by a linear mathematical model represented by f =S W V. 
5 A non- linear logical value S WL of 0 or 1 from the logical 
decision section is substituted for S w of this equation. 
The logical decision section receives an estimated 
observation quantity of the subsequent sampling time in 
form of a state quantity to be decided, and generates a 

10 logical value to decide the next model state. Substitution 
of the logical value for the switch element S w of the 
state selection section causes the output state quantity to 
offer f = 0 when the logical value is 0 , and f = v when the 
logical value is 1 . 

15 In the logical decision section of Fig. 46, U 

denotes an O R logical operator, H denotes an AND 
logical operator. Consequently, according to the logical 
expression of the logical decision section, when it is 
assumed that the logical 1 is given for the state 

20 quantities v x and f i being 0 or more, and the logical 0 
for others, and S x is logical values 0 or 1 , the logical 
value 1 is generated, when S x or f i is 1 and v x is 1, so 
that output state quantity f and input state quantity v of 
the state quantity selection section are connected with one 

25 another. 

(Explanation of example) 

8. Example of modeling for back-lash or bottom 
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projection 

Figure 47 is a view showing a typical structural 
model of a back -lash. Figure 47 shows an example in which 
a movable member A having a range of + L 0 is incorporated 
5 into an outer cylinder B. Energy by velocity v x and force 
f x is applied to the movable member A. Likely velocity V 2 
and force f 2 are applied to the outer cylinder B. 
Incidentally, y in Fig. 47 denotes a relative position of 
the movable member A wherein the center of the outer 

10 cylinder B is defined as 0 . 

According to the structural model in Fig. 47, at 
the range of ±L 0 for the movable member A, the movable 
member A and the outer cylinder B operate independently of 
one another, and in exceeding the range of +L 0 , both the 

15 movable member A and the outer cylinder B are united and 
operate as a rigid body. Consequently, it is considered 
that the back-lash (bottom projection) is associated with a 
basic function in connection and disconnection by a 
stiffness in which both the movable member A and the outer 

20 cylinder B are united. Rearrangement of this relation to 
establish modeling provides a functional model of Fig. 48. 
( 1 ) Functional model 

In Fig. 48, K g denotes stiffness when contacted; C 
g coefficient of viscosity resistance; v L velocity 
25 difference between the movable member A and the outer 

cylinder B; and x internal state quantity. In Fig. 48 (a), 
modeling is implemented from the mechanistic model through 
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a method of substitution of a non-linear stiffness. In Fig. 
48 (b), modeling is performed in form of a structural non- 
linearity in which a switch element is incorporated so that 
K g is connected or disconnected. It is possible for the 
back- lash, as mentioned above, to implement modeling in 
accordance with any method of the mechanistic model and the 
structural non-linearity. However, the basic function of 
the back- lash is to connect and disconnect stiffness of the 
outer cylinder A and the movable member B. Therefore, the 
modeling according to the structure non-linearity of part 
(b) of Fig. 48 is more desirably. The reason why it is so 
is that in case of representation of the mechanistic model 
of part (a) of Fig. 48, no switch element is incorporated 
into the functional model and the government equation; and 
thus it is difficult to identify the structural non- 
linearity. 

- - First, the functional model of the part (a) of Fig. 

48 is expressed by the following government equation. 



1 

' K g 
0 

0 
0 



0 



(7 1) 



0 1-1 

Next, the functional model of the part (b) of Fig. 
48 is expressed by the government equation and coupling 
condition equation as follows . 
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(7 2 ) 



(7 3) 



(2) Mechanistic model ( clearance model) 
The mechanistic model in part (a) of Fig. 48 is 
expressed by the following equation. 



L m =y- L o 

if dLj*0) then {K g =K g0 ) else (K g =0) 



(7 4) 



In eq. (74), the first line denotes a relative 
distance y of the movable member A; the second line a 
distance L ra between a wall of the external cylinder B and 
the movable member A; and the third line a stiffness K g 
of connection and disconnection. 

The mechanistic model of part (b) of Fig. 48 can 
be expressed by the following equation. 
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y =fv d _dt 

L m =y~ L o ( 75) 

if (\L m \ s> 0) then (V = 1) else (S^ = 0) 

In eq. (75), the first and second, lines are the 
same as eq. (74), and the third line denotes a condition 
decision equation of deciding a contact of the external 
5 cylinder B with the movable member A. In the condition 

decision, the contact state is S„y= 1 , and the non-contact 
state is S WY =0. Stiffness K g is connected or disconnected 
in accordance with an operation of the switch element S WY 
according to the condition decision Swy . 

10 (3) Result of simulation 

As one of examples, as shown in Fig. 49, let us 
consider a bottom projection of an impact absorption damper 
according to a coefficient C g of viscosity resistance as 
to a cylinder in which liquid is enclosed and a piston 

15 having a small aperture. In Fig. 49, M a denotes mass of the 
piston; C a a coefficient of viscosity resistance of the 
piston; M b mass of the cylinder ; and C b a coefficient of 
viscosity resistance of cylinder. Operational forces f a 
and f b of the piston and the cylinder are effected on the 

20 exterior to output velocities v a and v b . Regarding others, 
it is the same as the functional model shown in Fig. 49. 

Modeling of the structural model shown in Fig. 49 
can be implemented by adding M a * C a • M b • C b in Fig. 49 to 
the functional model of part (b) of Fig. 48 in form of the 
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additional function. Thus, the impact absorption damper is 
represented by a functional model shown in Fig. 50 wherein 
these additional functions are added to the functional 
model of part (b) of Fig. 48. Incidentally, x denotes 
internal state quantity. 

Coupling conditions of adding additional functions 
to both sides in Fig. 50 are expressed by the following 
equation. 



When the functional model and the cylinder of the 
piston and back-lash are united with eq. (76) and eq. (77), 
the following equation is obtained. 




(76) 




(7 7) 
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The functional and mechanistic models of an bottom 
projection damper shown in Fig. 50 were simulated with 
characteristic values shown in table 8. Figure 51 shows a 
simulation result. 

5 



10 



15 



20 
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In the simulation result of Fig. 51, velocity v a 
is reversed in positive and negative in order to compare 
with velocity v b . The simulation of Fig. 51 is to 
reproduce a state that the cylinder operates in accordance 
5 with a motion of the piston to which an operational force 
is applied. As the operation, 1 [N] was applied to the 
operational force f a of the piston, and after 0. 15 [s e 
c] the operational force —1 [N] in opposite direction 
was applied, and finally after 0. 3 [sec] the 

10 operational force was 0." At that time, the operational 

force f b [N] of the cylinder side was not applied. In the 
simulation result, from the top of Fig. 51 in order, 
f a [N] denotes an operational force which was applied to 
the piston; v a [m/ s e c] velocity of the piston; v b - Cm/ 

15 sec] velocity of the cylinder; f i * f 2 [N] force effected 
between the piston and the cylinder; and y [m] a relative 
position of the piston with respect to the cylinder. 
Incidentally, regarding the velocity on the second from the 
top, the solid line denotes the velocity v a of the piston 

20 at the driving side, and the dotted line denotes the 
velocity v b of the cylinder at the driven side. 

In Fig. 51, from the velocities v a • v b of the 
piston and the cylinder, they are in contact with each 
other about 0 . 0 3 [sec] and converges while the bottom 

25 projection is repeated, and finally after about 0. 0 8 [s 
e c] , they are united. It is understood that at the 
instance of the contact, a repulsion force generates by the 



140 



impact of the collision, and they converge while contact 
and separation are repeated. This is also understood from 
the variation of the impact load of force f i acting 
between the piston and the cylinder and the relative 
5 position y. Next, when the operational force f a is 

reversed, the piston contacts with a wall of the opposite 
side of the cylinder to show the same result. Finally, 
when the operational force f a is removed, the piston is 
united with the cylinder and stops . 

10 9. Friction clutch 

A clutch is a basic mechanical element for 
smoothly connecting and disconnecting two systems . 
Particularly, a friction clutch mechanism is well known as 
a mechanical clutch for changing over a system in operation. 

15 Here, modeling on the basic function of such a clutch will 
be considered. 

In a clutch, there are two types of clutch one of 
which is concerned with a normally connection scheme and 
another a normally disconnection scheme. Fig. 52 typically 

20 shows by way of example a structure of a clutch of a 

normally disconnection scheme in which a clutch piece of 
white squares at the right side is pressed to the left by 
press force F f of the clutch, so that a power transmission 
is performed between right and left systems. Applied to 

25 the clutch is rotational energy of (0 ± • T x from the system 
of the right side, and likely rotational energy of 0) 2 ' T 2 
from the right side. According to the clutch in which the 
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two systems are connected, the right and left systems are 
united with clutch stiffness K so that the same energy is 
transferred between both the systems. In the sliding state, 
the right and left systems are connected with one another 
5 by the frictional torque of the clutch so that an angular 
velocity difference occurs between the right and left 
systems . Loss caused by the angular velocity difference 
and the frictional torque is consumed in form of thermal 
energy. In the state in which the clutch is disconnected, 
10 the right and left systems are separated and independent of 
one another . 

(1) Basic function 

Next, the basic function of a clutch will be 
considered. The function of a clutch can be divided to the 
15 original functions of connection, sliding and disconnecting, 
and the additional functions of parameters possessed by two 
systems . The original functions divided to the three 
states can be considered as follows. 

® Connection: A state that no difference in 
20 velocity exists between two systems so that the clutch does 
not slip. 

(2) Sliding: A state that two systems are driven by 
a frictional torque of the clutch. 

(D Disconnection: A state that the frictional 
25 torque is 0 so that two systems are separated. 

From these three states, it is understood that 
sliding by the friction torque controls connection and 
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disconnection of the clutch. This relation is shown in Fig. 
53 . 

Figure 53 is a view showing torque transfer 
characteristics of a clutch, where the horizontal axis 
5 denotes friction torque T f and the vertical axis denotes 
load torque T x • T 2 . T f0 shown in Fig. 53 denotes a 
transfer torque of the clutch by friction. The transfer 
torque T f0 causes the right and left systems to be united 
with clutch stiffness K in the area wherein the load torque 

10 T x • T 2 is small, so that the clutch offers the connecting 
state. In other area, the load torque Tj. • T 2 is supplied 
by the transfer torque T f0 , so that the clutch offers the 
sliding state. The transfer torque T f0 varies in 
accordance with a compression force F f of the clutch, 

15 wherein when T f0 = 0 , load torque T 1 =T 2 =0, so that the 
clutch offers the disconnection state. Therefore, a 
relation between the connection and the sliding in Fig. 53 
is controlled in accordance with a relation in magnitude 
between the load torque T x * T 2 to be applied via the 

20 clutch to two systems and the friction torque T f . For 

example, in the event that the load torque T 1 • T 2 is less 
than friction torque T f0 , the connection is established, 
and in the event that the load torque T : • T 2 is more than 
friction torque T f0 , the sliding is generated owing to the 

25 friction torque T fo . From this operation, the clutch is 
concerned with the basic function such as the clutch 
stiffness K with which two systems are united, the 
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transfer torque T fo controlling the sliding of the clutch, 
and change over of both the systems. 

(2) Functional model (clutch model) 
Fig. 54 shows functional and mechanistic models of 
5 a clutch which are subjected to modeling in accordance with 
the basic function of the clutch. Fig. 54 shows a clutch 
model of a scheme of normally disconnecting. Part (a) of 
Fig. 54 shows a clutch model connecting by the clutch 
stiffness K. Part (b) of Fig. 54 shows a clutch model in 

10 which the clutch stiffness K is replaced by the 

coefficient D of viscosity resistance to perform a low 
leveling. In Fig. 54, K denotes stiffness of the clutch 
at the time of connection; D a coefficient of viscosity 
resistance of a frictional piece of the clutch; T f friction 

15 torque; T c _ stiffness torque of the clutch; and S wc a 
switch element for selecting between T c _ and T f . A 
mechanistic model for operating the clutch is incorporated 
into a functional model constituted of those elements, x 
denotes an internal state quantity. The clutch model of 

20 part (b) of Fig. 54, in which the clutch stiffness K is 
replaced by the coefficient D of viscosity resistance, 
generates a slight slip in the connecting state. While the 
general characteristic of the coefficient D of viscosity 
resistance is of linear, transforming the characteristic 

25 into a non-linearity by the mechanistic model to provide a 
negative coefficient of viscosity resistance makes it 
possible to reproduce phenomena such as abnormal vibration 
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and abnormal sound (squeak). 

Government equation of part (a) of Fig. 54 is 
expressed by the following equation. Incidentally, the 
connecting condition of the mechanistic model is omitted. 
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In eq. (79), the first line denotes a state 
equation; second and third input and output equations; and 
fourth an observation equation of the clutch torque. 
Likely, the government equation of part (b) of Fig. 54 is 
expressed by the following equation in which the state 
equation and the internal state quantity are removed from 
eq. (79). 
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(3) Mechanistic model (coupling decision model) 
Parts (a) and (b) of Fig. 54 are substantially the 
same mechanistic model which outputs condition decision S wc 
wherein stiffness torque T c _ of the clutch is given in 
form of an estimated observation quantity, and transfer 
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torque T f0 in form of substitution quantity. The contents 
of the mechanistic model is as follows . 

The transfer torque T f0 can be determined from 
equation of the following friction torque. 



-^f A f R f F flh (81) 



Where M f denotes a coefficient of kinetic friction 
of a friction material; A f a sliding area of a frictional 
surface; R f an average radius of the frictional surface; 
and F f a clutch compression load, and the positive and 
negative of T f0 is coincident with an estimated 
observation quantity T c _. 

According to this example , there is shown a model 
in which a friction torque T £ of the clutch is represented 
by a kinetic friction torque. Consequently, when the 
behavior of the clutch is considered in details, it is 
necessary to consider a variation of a friction torque due 
to the static friction, the kinetic friction and the 
frictional heat. 

Next, the condition decision for deciding a 
sliding of the clutch and the connection state is expressed 
by the following equation. 

iy(|l) 0 |*|r c |) then (S wc = 1) else (S wc =0) 

(8 2) 
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In eq. (82), when transfer torque T f0 is not less 
than stiffness torque T c , S wc = 1 is set up to establish a 
connecting state, and in opposite S wc =0 a sliding state. 
In the sliding state, an internal state quantity X k of the 
5 dynamic storage element is initialized in accordance with 
the following equation. 

if (S wc = 0) then (x k = 0) (83) 

According to the clutch model of a normally 
10 connecting scheme, there is given an equation 

I T f0 I ^ I T c _ I in which a decision of comparison operators 
of the condition decision equation (82) is reversed, and a 
switch element S wc 0 is interchanged by S wc _i . 
(4) Result of simulation 
15 Drive systems and load systems of moments of 

inertia J D * J l and coefficients of viscosity resistance 
*C D • C L are added to right and left of the structural model 
shown in Fig. 52, and this is connected to the exterior 
with the input torque T D • T L and the state quantity of the 
20 output angular velocity co D • co L . The functional model is 
shown in Fig. 55. 

For the functional model of Fig. 55, there is used 
a clutch model taking into consideration the clutch 
stiffness of part (a) of Fig. 54. x denotes the internal 
25 state quantity. The government equation of Fig. 55 is 
expressed by the following equation. 
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Fig. 56 shows a result of a simulation in which 
characteristic values of table 9 are used for the 
functional model of Fig. 55. Incidentally, in Fig. 56, 
angular velocity co D is reversed in positive and negative 
for the purpose of comparison with the velocity co L . 
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According to the simulation of Fig. 56, 0. 1 [N] 
of torque T D is applied to a drive system, and a clutch 
friction torque T f is raised at 0. 2 5 [Nra/ sec] to 
connect the clutch. After the connection, after 1. 5 [s e 
5 C ] , load torque T L = 0 . 14 5 [Nm] is applied during 

0. 5 [sec] on a step basis to slide the clutch, and at 
the final 3. 5 [sec] the clutch is released. In Fig. 
56, from the top in order, T D [Nm] denotes torque of a 
driving system; T f [Nm] (a side load) friction torque of 

10 the clutch; T L [Nm] torque of a load system; co D drive 

angular velocity; co L load angular velocity; and T L torque 
of the clutch. 

In Fig. 56, in activation, the angular velocity 0) D 
of the drive system rises to increase the friction torque 

15 T f so that the load angular velocity co L of the load 

system rises too. However, the driving system is pulled in 
the load system halfway so that the angular velocity co D 
slow downs, and both the driving system and the load system 
are coincident with one another at the point of about 0 . 5 

20 [ S e C ] to offer a connecting state (engagement), so that 

they are united to rise in rotation. Up to the coupling, 
the load system is driven by the friction torque T f , and 
the transfer torque T ± of the clutch and the friction 
torque T f are the same as one another. After the 

25 engagement, the transfer torque T 1 of the clutch is 

replaced by torque by the clutch stiffness, and rapidly 
lowered . 



151 



Next, when torque T L is applied from the exterior 
to the load system, the torque applied to the clutch 
exceeds the friction torque T f so that the clutch begins 
to slide whereby an angular velocity difference is 
5 generated between co D and co L . The load torque T x of the 
clutch at that time is again replaced by the friction 
torque T f . Removal of the load torque T x causes the 
angular velocity difference between C0 D and CO L to disappear, 
so that the clutch connecting state is established whereby 

10 the rotation rises again. Finally, when the friction 
torque T f is set to 0 to disconnect the clutch, the 
driving system and the load system are individually 
decelerated and stopped. 

10. Friction brake 

15 A brake is one of the basic mechanical elements 

for a velocity control by compulsory absorption of energy, 
and for absorbing a storage energy such as moment of 
inertia. There are many types of brake. Here a brake 
utilizing a frictional force will be considered. 

20 A friction brake controls, in a similar fashion to 

that of the clutch, a braking torque by a frictional force. 
Application of the clutch model to a brake can be 
implemented in such a manner that the rotating driving 
system are fixed on one side of the load system so that an 

25 energy consumed in sliding of the clutch is supplied to the 
driving system in form of a braking energy. Regarding a 
fixing at the load side, there are two ways. One of the 
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two ways is a method in which stiffness is connected to the 
load side, or alternatively moment M L of inertia of Fig. 55 
is enlarged. This is applied when an influence of the 
braking force on a mounting section is considered too. 
5 Another is a method in which moment M L of inertia and 

coefficient C L of viscosity resistance at the load side of 
Fig. 55 are removed and the angular velocity 0) 2 is set to 
0 . This method is applied when the braking function of the 
brake is considered. Here, the later brake model is 

10 considered. In view of the foregoing, modeling of the 

brake is effected hereinafter. Fig. 57 shows typically a 
structure model of a brake. In Fig. 57, the angular 
velocity co D and drive torque T D are added to the drive 
side of the shaft, and the angular velocity co D and drive 

15 torque T D are applied to the load side. Friction torque 
T F of the brake is effected in form of the braking torque 
on the shaft. 

(1) Functional model (brake model) 

A brake is different from a clutch and has no 

20 function of transferring motive power, and thus serves as a 
model for applying a braking energy in form of a load to a 
braking shaft to transfer the motive power. Consequently, 
the brake model is a functional model in which a clutch 
model is placed in parallel with the model of the drive 

25 shaft. Fig. 58 shows a functional and mechanistic model of 
the brake . 

In Fig. 58, moment of inertia J D and coefficient 
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of viscosity resistance C D are added to the driving side of 
the driving shaft, and moment of inertia J L and 
coefficient of viscosity C L resistance are added to the 
load side. In Fig. 58, «i and T 1 denote angular velocity 
5 and braking torque. With respect to function of generating 
the braking energy, input and output state quantities <x> 2 
and T 2 of the right side of Fig. 54 are removed from part 
(a) of Fig. 54. While the mechanistic model is almost the 
same as part (a) of Fig. 54 in characteristics, £t f is 

10 replaced by a coefficient of kinetic friction of the brake 
pad; A f a sliding area; R f a sliding radius; F f a brake 
operational load. 

Inertia moments J D and J l at the drive side and 
the load side, which are shown in Fig. 58, are the same as 

15 internal states x d ' and X L , respectively. Thus, it is 
possible to derive inertia moments J D + J l through 
addition and connection of the differential quantity. The 
government equation, in which inertia moments J D and J L 
are combined, is expressed by the following equation. 



25 
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In eq. (85), the first and second lines denote 
state equations; the third line an input and output 
equation; the fourth and fifth lines observation equations 
and the sixth line an estimated observation equation. 
5 Regarding the mechanistic model, it is the same as eqs . 
(81) to (83) of the clutch model. 

(2) Result of simulation 

Fig. 59 shows a result of a simulation in which 
characteristic values in table 10 are applied to the 
10 functional and mechanistic models shown in Fig. 58. 
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In Fig. 59, from the top in order, co D [rad/se 
c] denotes an angular velocity of the driving shaft; T f 
[Nm] a control torque; T c _ [Nm] an estimated observation 
quantity; T D [Nm] driving torque; and S 

wc condition 

5 decision. 

In execution of the simulation, at the same time 
of the actuation, driving torque T D = 1 [Nm] is applied, 
and at 0. 0 2 [sec]a brake operational load F f =8Xl 0 4 
[N] is applied to effect a brake, and at 0. 0 6 [sec] 

10 the torque is changed over to a reverse driving torque T D = 
— 1 [Nm] . And at 0. 0 8 [sec] the brake operational 
load slow downs at the velocity of F £ =2 X 1 0 4 [N/s e c] 
and at 0. 1 [sec] the brake is disenabled. 

In Fig. 59, when the brake is effected at 0. "0 2 

15 [sec], the state of T f <T c _ appears and thus the 

friction torque T f causes the brake torque to effect so 
tbrat an angular velocity CO D is decelerated. In the state 
of T f >Tc_, it is stopped compulsively by the brake. When 
the driving torque is reversed at 0. 06 [sec] , the 

20 reaction causes a slip of the brake to instantaneously 
appear on the angular velocity co D so that condition 
decision S wc immediately decides. When the brake is 
weakened gradually from the time point of 0 . 0 8 [sec] , 
the brake begins to slip from the time point of T f <T c _ , 

2 5 and the angular velocity co D is also increased (the 

opposite direction). When the brake is disenabled, it is 
accelerated. 
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11. Automatic-reset mechanism 

An automatic-reset mechanism is, as shown in Fig. 
60, a mechanism wherein a return torque is applied, to part 
of a lever to return to a stop position by stiffness K R of 
5 a return spring. This lever is controlled in an operation 
range in accordance with stiffness K c and viscosity 
resistance D c of a stopper made of rubber. Functional and 
mechanistic models of this mechanism are shown in Fig. 61. 
In Fig. 61, a pair of angular velocity CO s and 
10 torque T s of a lever is connected to the exterior in form 
of input and output state quantity. With respect to the 
internal characteristic, the return spring is represented 
by stiffness K R , and the rubber-made stopper is 
represented by stiffness K c and coefficient D c of non- 
15 linear viscosity. Distances between the rotating center of 
the lever acting on this and the load points are 
represented by sizes L R and L c shown in Fig. 60, which 
are factors of Fig. 61. In Fig. 61, X denotes internal 
state quantity. 
20 First, the government equation of Fig. 61 is 

expressed by the following equation. 
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Equation (86) is a state equation of one raw-one 
column, wherein since springs of K R and K c of Fig. 61 are 
connected in parallel, they are united into stiffness of a 
single spring to perform a low leveling. 
5 Next, a mathematical model is determined on the 

mechanistic model for implementing the function of Fig. 61. 

An estimated rotating angle 8 R _ is controlled in 
movement in accordance with the lever operating range 6> max 
[r ad] shown in Fig. 60. This condition decision S wc is 
10 expressed by the following equation. 
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Regarding the condition decision of eq. (87), when 
condition of S wc = 1 is applied, the lever is compulsively 
stopped, and when condition of S wc =0 is applied, the 
lever rotates freely. 
5 Finally, it is understood through a trial 

examination that effecting force of the rubber-made stopper 
for compulsively stopping the lever is different between a 
compression direction and an expansion direction. So, this 
is represented by a non-linearity of a coefficient D c of 

10 damping resistance. It is estimated that this non- 
linearity is caused by the fact that while a load of the 
rubber effects in the compression direction, it does not so 
effect in the expansion direction because of a delay of the 
expansion. In fact, a compression stiffness K c of the 

15 rubber is also non-linear according to the deformation 

quantity, but this stiffness is regarded as linear modeling. 
Figure 62 shows non-linear characteristics of an 
attenuation resistance factor D c . 

In Fig. 62, a damping resistance coefficient in 

20 the compression direction in a state that the lever 

contacts to the stopper is denoted by D C H , and a damping 
resistance coefficient in the expansion direction is 
denoted byD CL , and it is assumed that the damping 
resistance coefficient is not varied by a quantity of 

25 deformation. Modeling is implemented in such a manner that 
the directional property of the damping resistance 
coefficient is decided through an association of the 
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forward rotation and backward rotation of the lever with 
the compression and expansion direction, and the damping 
resistance coefficients are changed over. The condition 
decision S WCD is expressed by the following equation. 
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In eq. (88), CO R _ determines a rotational direction 
of the forward rotation and the backward rotation, and when 
its product with 8 R _ is positive, the judge is compression, 
and when it is negative, the judge is expansion. Selection 
5 of the damping resistance coefficients D C _ H or D C _ L is 

performed by selection switch elements S WCD _ 0 and S^^! 
which simultaneously have no states of 1 and 0 . When the 
judge is not applied, the compression side is denoted by S 
WCD _o- 1 , and when the judge is applied, the expansion 
10 side is denoted by S WCD _ 1 = 1 . 

12. Rattle device 

A rattle device is a mechanism for preventing an 
interference between a manual operation from the exterior 
and an automatic operation by an actuator. Fig. 63 is a 

15 view showing a construction of a rattle device. 

In Fig. 63, a lever driven by a motor penetrates 
ari aperture having width L H provided on a slider and 
couples with the slider. In Fig. 63, X L denotes a moving 
quantity of the lever based on the lever stop position 

20 shown in Fig. 60; X s an estimated moving quantity of the 
slider; and Y F a relative distance between the lever and 
the slider. Where X L and X s are an estimated moving 
quantity. Incidentally, an operation of the rattle device 
will be described later. Fig. 64 shows the functional 

25 model of the rattle device in Fig. 63. 

In Fig. 64, K F denotes a compression stiffness of 
a slider pressing side; M s mass of a slider; D s a 
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coefficient of viscosity resistance; and L F a contact 
position of the slider with the lever shown in Fig. 60. 
Variation of L F due to the rotational angle of the lever 
is neglected. In Fig. 64, an angular velocity CO F and a 
5 torque T F are coupled with a decelerator mechanism, and an 
output velocity v s and a load F s are connected with an 
external door latch mechanism. X s denotes an internal 
state quantity. 

First, a government equation of Fig. 64 is 
10 expressed by the following equation. 
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The upper part of eq. (89) shows a government 
equation wherein the first and the second lines denote 
state equations, the third and the fourth lines denote 
input and output equations respectively, and the third line 
5 denotes an estimated observation. In eq. (89), the lower 
part shows an equation of an estimated storage difference 
quantity of mass M s wherein the first line denotes an 
internal storage observation quantity X s s , and the second 
line denotes an estimated moving quantity V s s of a slider. 
10 Incidentally, S amp _ t is a sampling period. 

A mathematical model of the functional model 
appearing at the upper side of Fig. 64 is introduced as 
follows . 

First, an estimated observation quantity V FF _ of a 
15 lever position is derived from a rotational angle of a 
deceleration function in accordance with the following 
equation. 

V FF_ = L F 6 F (9o> 

20 Next, with respect to a moving quantity X s of a 

slider, since a difference system of mass M s includes a 
moving distance V ss , the moving quantity X s is expressed 
by X S = V SS . A moving quantity X L of a lever is expressed 
by eq. (90). From the moving distance V ss and the eq. 

25 (90), the relative position between the lever and the 
slider is expressed by the following equation. 
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Y F ~ X S ~ X L (91, 

Fig. 63 shows the rattle device in the state that 
the slider is in a state of an unlock side, and the lever 
is pushed to a lock side. In Fig. 63, a view of the top 
shows the state that the lever starts to be pushed, a view 
of the middle shows the state that the lever stop position 
is located at the center of width L H of an aperture, and a 
view of the bottom shows the state that the lever is 
completely pushed to the lock side. From Fig. 63, it is 
understood that when the absolute value I Y F I of eq. (91) 
is within the range L H , the lever rattles and thus the 
slider does not move. The condition decision S WF of the 
rattle is expressed by the following equation. 

if (|y F | s ^L) then S m =1 else 5^=0 

(92) 

Finally, from eq. (92), on the condition 
establishment of S WF = 1 , the lever pushes the slider, and 
on condition failure of S WF =0, the lever rattles. The 
moving quantity X s of the slider at that time determines 
directions of push and rattle in accordance with an 
estimated velocity v s _ of the slider and the relative 
position Y F . The condition decision S 

wp expressed by 

the following equation. 



170 



if (y s Y F > 0) then = 1 else 

(9 3) 

Fig. 65 is a block diagram of a fourth non-linear 
5 characteristic reproducing apparatus according to an 

embodiment of the present invention. Fig. 66 is a view 
useful for understanding a determination state of the non- 
linear characteristic reproducing apparatus shown in Fig. 
65. In both the figures, regarding variables, ( k — 1 ) 

10 denotes the previous sampling time; (k) the present 

sampling time; and (k+1) the subsequent sampling time. 
+n represents n-order subsequent sampling time from the 
present time to the future. State quantity v taking k as 
the suffix represent is representative of observation 

15 quantity of the respective time. In the figures, Av 
denotes a variation width of the estimated observation 
quantity between the present sampling time and the 
subsequent sampling time; V s deviation between the present 
observation quantity v (k) and the decision value V TH ; t smp 

20 a sampling period. In Fig. 65, a stable state decision 

section provides such a decision that when deviation V s is 
smaller than the variation width A v , the decision value 
is given with 1, and when deviation V s is larger than the 
variation width A v , the decision value is given with 0 . 

25 According to the prior art, the decision is made 

through a matter that observation quantity is within a 
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range of the designated stable width (dead zone). For this 
reason, variation of sampling period and observation 
quantity notable in variation causes a stability decision 
to be unstable. To the contrary, according to the decision 
5 method proposed here, wherein no stable width (dead zone) 

is provided, it is possible to prevent non- stability due to 
variations of the sampling period and the state quantity. - 

As shown in Fig. 67, an adoption of a decision 
value V TH1 for rising and a decision value V 

th2 for going 

10 down may establish a hysteresis for rising and going down 

of the state quantity, so that the more stable decision can 
be expected. 

Fig. 68 is a principle explanatory view of a fifth 
non-linear characteristic reproducing apparatus according 

15 to an embodiment of the present invention. Incidentally, 
Fig. 68 is a view useful for easily understanding a 
principle the fifth non- linear characteristic reproducing 
apparatus according to an embodiment of the present 
invention, and the fifth non-linear characteristic 

20 reproducing apparatus according to an embodiment of the 
present invention is not confined to Fig. 68 and the 
associated description. 

In Fig. 68, a stationary state of variation value, 
which is generated in accordance with state quantities S 2 

25 and S 3 , is fed to a normalization response section having 
a response range of 0 to 1 in which the stationary state of 
variation value is reproduced in form of a non -linear 
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response value varying on a time basis in accordance with a 
non-linear response characteristic of the normalization 
response section. A result of reproduction is fed to a 
characteristic generation section in which it is 
5 transformed to a variation of characteristic R of a state 
quantity transformation section and substituted, so that a 
slow response characteristic is reproduced (reflected) in 
the whole system of linear model. 

The use of the apparatus for reproducing the time 

10 history of response makes it possible to remove an inherent 
value dominating an influence of the temperature rise due 
to the energy loss concerned with the slow behavior and an 
influence for variations of the environment such as the 
outside air temperature and the atmospheric pressure, from 

15 an inherent value dominating a complicated transient 

characteristic of the linear model of the whole system. 
Here, the range of the removable slow transient response is 
restricted to the range in which inherent values of an 
electric system and a mechanical system are not stimulated. 

20 (Explanation of example) 

13. Example of modeling for temperature rise 
Here, let us try to establish modeling in which a 
relatively slow transient change is reproduced in form of a 
non-linearity. 

25 (1) Normalized primary response model 

In many case, it happens that many of 
characteristics varying with the passage of time is 
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expressed in form of primary or secondary response 
characteristic. For this reason, there is used a discrete 
equation in which the primary and the secondary response 
delays are normalized in a range of 0 to 1 and are 
considered as factors for non-unit quantities. The general 
equation of the primary response is expressed by the 
following equation. 

z uk + i) = Pi z (k) + 1i u i(k-> 

=CiZ(k)+ d i u i 

In eq. (94), P x denotes an inherent of a discrete 
primary delay, which converges in a range of 0 1 and 
disperses in a range of more than 1 . And r 1 denotes a 
time constant of the system; t smp a sampling period on a 
discrete basis; and suffix k a sampling time. 

In eq. (94), z 1(k+1) and z 1(k) denote internal state 
quantities; U 1(k) and Yim input and output state 
quantities; Qi a coefficient of input state quantity; and 
d a coefficient of output state quantity. In the 
stationary state of eq. (94), Z 1(k+1) = Z 1(k) = U 1(k) are given, 
and thus there is generated no difference between the input 
state quantity and the output state quantity. Therefore, 



(9 4) 
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the state equation of the upper side of eq. (94) can be 
normalized when the coef f icient q x of input state quantity 
is expressed by the following equation. 

5 <?i = (l-^i) (9 6) 

The primary response (normalized) of the 
subsequent sampling period, wherein eq- (96) is substituted 
for Qi of eq. (94), the coefficients of the output 
equation are given with C x = 1 • d x = 0 , and the output 
10 equation is omitted, is expressed by the following equation. 

Z l(£+1) = Pl Z l{k) + (1 ~ Pl) U l (97) 

(2) Normalized secondary response model 
15 Likely, when the secondary response is regarded 

also as an interference of two primary responses and is 
rearranged, the secondary response is expressed by the 
following equation. Incidentally, the rearranging process 
is omitted. 

20 

Z a {k+ 1) = Pa Z a (k) +Q-~Pa )( X " 8 ab 

Z b {k+ i ) = P h z Hk) + (1 - P b )(l + g ab )u m \ 

i \ (98) 

Where p a and p b denote inherent values associated 
with eq. (95), and g a b denotes a coefficient indicative of 
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a degree of the interference between two primary responses . 
They are expressed by the following equation. 



Pa = expH^ /rj 
Pb =exp{-^ /rj 

'8ab = ( T a +T J/( T « ~ X b) 



(9 9) 



In eq. (99), r a and r b denote tie constant, and 
t smp denotes a sampling period. 

(3) Temperature rising model 

By way of example, a temperature dependency of a 
winding resistance of a motor is subjected to a modeling. 
Application of an electric power to a motor causes a 
temperature to rise owing to self -heating by the consumed 
power P M • The temperature rise T 0 is expressed by the 
following equation, where 6 M is a coefficient of heat 
radiation of the motor. 



determined by the consumed power of coil resistance P M - 
As the temperature rises and the current is reduced, the 
consumed power is reduced. And as a result, the mechanical 
output is reduced. When the motor current is denoted by I 
M and the coil resistance is denoted by Rm, the consumed 




M 



(10 0) 



From eq. (100) the temperature rise T v is 
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power P M is expressed bt the following equation. 



(10 1) 



It is known that a process wherein the temperature 
generated in eq. (100) becomes a stationary state offers a 
primary delay of response. Thus, the response of the 
temperature rise wherein the temperature rise of eq. (100) 
is set up to the target value is expressed by the following 
equation in accordance with eqs . (97) and (100). Where r M 
in the equation denotes a thermal time constant. 



M 



(102) 



Next, a coil resistance R M in coil temperature T x 
having characteristics of a temperature coefficient a w , 
relation between a reference temperature T MS and a 
reference resistance R MS is expressed by the following 
equation. Where T x denotes the estimated temperature T 
x(K +i) represented by eq. (102). 

R M = R MS ^ + &W ( T x(k+1) ~ T MS ) j 

(103) 

Equation (103) gives a coil resistance value of 
the subsequent sampling time. When a atmospheric 
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temperature is considered, the atmospheric temperature is 
added in the term of the temperature of eq. (103). 
(4) Functional model 

Fig. 69 shows a functional model incorporating 
thereinto a mechanistic model in which the above-mentioned 
relation is rearranged to establish a modeling. 

The functional model of Fig. 69 is a combination 
of a supply element and a basic functional element of a 
fluid loss properties . The functional model supplies the 
consumed power P M of a coil resistance of a motor to the 
mechanistic model and feeds from the mechanistic model an 
alteration value of a coil resistance R M due to the 
temperature rise to the functional model. The output 
difference quantity of the coil resistance R M stores the 
temperature rise T x in form of a temperature energy. 

The functional model of Fig. 69 is expressed by 
the following government equation. 



' R B R b 
1 E n 



(104) 



(5) Result of simulation 

Fig. 70 shows a result of simulation on the 
functional model of Fig. 69 in accordance with the 
characteristics of table 11. 
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In the simulation, the voltage of 12 [V] is 
applied at 0 [sec], and at 1500 [sec] it is turned off. At 
the later time, assuming a lock state in which a rotation 
of the motor is fixed, an input voltage V M1 is set to 0 [V] . 
5 The reason why in Fig. 70, the target value T 0 is lowered 
owing to the temperature rise after the power source is 
turned on is that the temperature rise involves the 
increment of the resistance value of the coil and the 
consumed power is lowered. With the passage of time T v 
10 and T x are coincident with each other and it offers a 

stationary state wherein heating and heat radiation are 
balanced. 

Table 11 

15 Characteristic table of coil temperature rise 



Characteristic name 


Symbol 


unit 


Characteristic 
value 


Coil resistance of reference 
temperature 

Reference temperature of coil 
resistance 

Coefficient of coil temperature 


R MS 

a w 


[Q] 

C°c] 
[sec] 


2.0 
25.0 

4.3X10-* 


Thermal time constant of motor 
coil 

Coefficient of heat radiation 
of motor coil 


t M 


CW/deg] 


200.2 
1.8 


Source voltage 


E 0 


[V] 


12.0 


Internal resistance of source 


Rb 


[□] 


0.001 



14. Positive characteristic thermistor 

(1) Outline of overload protection function 

A motor is subjected to an overload and a 
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compulsory stop by an equipment connected to an output of a 
mechanical system. At that time an overcurrent nearly as 
large as the starting current conducts through a coil 
resistance of the motor. In the event that this state 
5 continues a long time, the motor is heated and burned. 

Regarding the overload, since torque and current of the 
motor is in proportion to a motor constant M M , a 
provision of the function for detecting and suppressing 
overcurrent makes it possible to prevent such a burning of 

10 the motor. As methods for protection of the motor from 
burning, there are many ways such as the simple use of a 
fuse fused by an overcurrent and a complicated way on an 
electronic basis. Here, there will be considered a method 
of preventing burning in such a way that a single positive 

15 characteristic resistance device, which is usually used in 
a miniature motor for a car, is adopted to detect the 
overcurrent. Here, the positive characteristic resistance 
device is referred to as a positive characteristic 
thermistor. Fig. 71 is a perspective view of the positive 

20 characteristic thermistor. 

The positive characteristic thermistor has such a 
non-linear characteristic that when a temperature is low, a 
resistance value is small, and when the temperature is high, 
the resistance value is large. According to such a 

25 characteristic, when the motor current is within a range of 
tolerance, the resistance value of the positive 
characteristic thermistor is small, and when the motor is 
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in an overload state, the resistance value is increased 
owing to self -heating due to the overcurrent so that the 
current is suppressed. This working is implemented by a 
contactless switch in which the non-linear resistance 
5 characteristic of the positive characteristic thermistor 
having temperature dependency properties is utilized to 
perform detection and blocking of the overcurrent . 

(2) Functional model of an overload protection 
Fig. 72 is a view showing a functional model of an 

10 overload protection in which the positive characteristic 
thermistor is connected in series to the motor so that 
motor current I M ( I x ) directly conducts. 

In Fig. 72, a resistance R P denotes a non-linear 
resistance of a positive characteristic thermistor; t P a 

15 temperature of a static characteristic thermistor; and t F 
an ambient temperature. According to the mechanistic model 
of the overload protection characteristic, there is 
generated a non-linear resistance in which the positive 
characteristic thermistor detects an overcurrent to be 

20 suppressed, and the non-linear resistance thus generated is 
substituted for R P . Further, according to this 
mechanistic model, V T and I T are connected to the motor, 
and V P and I P are connected to the operational switch at 
the power supply side. 

25 The functional model of the overload protection 

shown in Fig. 72 is expressed by the following government 
equation . 



181 



■R P -1 





V 







(10 5) 



(3) Mechanistic model of a positive characteristic 
resistance element 

A non-linear resistance characteristic of a 
positive characteristic thermistor for performing an 
overload protection of the motor can be represented by 
temperature - resistance ratio characteristics shown in 
Fig. 73. 

In Fig. 73, with respect to the relation between 
the temperature t P of the positive characteristic 
thermistor and the resistance ratio, the resistance ratio 
is represented by logarithms, since variations at the high 
temperature side are remarkable as compared with the low 
temperature side. The characteristic is defined in such a 
manner that the resistance ratio of the reference 
temperature of 25 [X] is given as 1, and a temperature, 
wherein the resistance ratio becomes twice, is given as a 
resistance critical point. In Fig. 73, an area shown with 
a hatching portion denotes a tolerance of the positive 
characteristic thermistor to be used. At that time, the 
temperature characteristic of the positive characteristic 
thermistor offers, as shown in Fig. 73, such a compound 
characteristic that at the low temperature side a negative 
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characteristic appears and at the high temperature above 
the resistance critical point or so a positive 
characteristic appears- To establish modeling, the 
modeling is expressed by an approximation equation in which 
5 the compounded two characteristics are combined. 

The resistance ratio at the low temperature side 
is expressed by the following equation from the general 
equation of the negative characteristic thermistor 
represented by the inverse of the temperature. 

10 

^ =exp { B -(^h;-^h:)} < 106 > 

The resistance ratio at the high temperature side 
is expressed by the following equation which is an 
15 approximation equation of the positive characteristic 

wherein the temperature term shown in eq. (106) is replaced 
by a proportion from an inverse proportion. 

^-" =exp t^ _7W) } 

20 From eqs . (106) and (107), the approximation 

equation of the non-linear resistance of the positive 
characteristic thermistor is expressed by the following 
equation . 
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R P =R PS {ip P _ L +ip P _ H ) (108) 

where t P : a temperature ['C] of the positive 
characteristic thermistor 

R P : a resistance value [Q] at a temperature t x 
5 R PS : a reference resistance value at the 

reference temperature 2 5 [°C] 

"5 r P _ I ^ f P _ H resistance ratios at the low 
temperature side and the high temperature side, 
respectively 

10 T PL a reference temperature [°C] in which the 

resistance ratio is 1 

T PH : a temperature [°C] at the critical point of 
resistance of a transition to the high temperature side 
B PL : a constant [K] corresponding to B- 
15 characteristic of the negative characteristic thermistor 
B PH : a constant [K] corresponding to B- 
charact eristic of the positive characteristic thermistor 

Table 12 shows characteristic values of the 
positive characteristic thermistor which some device maker 
20 opens to the public. The resistance characteristic applied 
to eqs. (106) and (107) is indicated with the hatching 
heavy line in Fig. 73. 
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In Fig. 73, a resistance variation is small up to 
100 [°C] or so of the resistance critical point in which the 
motor coil is not damaged, and a low resistance appears. 
When the overload current conducts through the motor, a 
5 consumed power P P of the positive characteristic 

thermistor resistance is increased, so that a self -heating 
occurs owing to Joule heat. The self -heating involves 
further temperature rise t PU of the positive 
characteristic thermistor, so that the resistance value 

10 rapidly increases at the resistance critical point or so 
and whereby the heating further advances . On the other 
hand, the rapid increment of the resistance value serves to 
decrease the motor current , so that the temperature 
stabilizes at the balanced point of heating and heat 

15 radiation of the positive characteristic thermistor and 
settles down at a temperature t P . At that time, the 
current maintains a smaller value with respect to the upper 
limit in which the motor is subjected to burning. 

With respect to the self-heating and heat 

20 radiation of the positive characteristic thermistor which 
dominates the overload protection, in a similar fashion to 
that of the above-mentioned motor coil, the temperature 
rise t PU [deg] is determined the heat radiation constant 
(5 P [W/deg] . The temperature rise t PU [deg] still in 

25 progress is determined in accordance with the thermal 

constant x p of the positive characteristic thermistor. One 
in which the ambient temperature t F [°C] is added to the 
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temperature rise t PU [deg] is a temperature t P [°C] to be 
stored in the positive characteristic thermistor. t F 
denotes the atmosphere temperature [°C] . Fig. 74 shows a 
relation of those in form of the mechanistic model. 
5 (4) Stationary functional model of a motor-lock 

state 

Now let us consider an overcurrent of the motor- 
lock state and the behavior of the positive characteristic 
thermistor in the event that a rotation of the motor is 

10 compulsively stopped, and verify the utility of the 
overload protection characteristic. Fig. 75 shows a 
stationary functional model in which the condition of the 
motor- lock state is incorporated. 

In Fig. 75, R P • R M are varied in resistance value 

15 by self -heating by a lock current I 0 _ L , but not varied in 
another. In Fig. 75, with respect to the brush potential 
drop E MB , since the motor is in the stop state, E MB = 0 is 
to be set up. However, in view of the fact that even in 
the stop state the armature is associated with the micro 

20 vibration, the brush potential drop is considered. Further, 
in Fig. 75, since the induced voltage V w is not generated, 
electric resistance elements, which are connected in series 
with each other, are directly connected to the battery, so 
that very large current I 0 _ L conducts through the motor. 

25 The mathematical model of the functional model of 

Fig. 7 5 is expressed by the following equation. 
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j Bp E MB 

°~ L R + R + R d0 9) 



In eq. (109), since the power supply is in the 
state of turn-on, the operational switches set up as S WE 

= S WA = 1 • 

(5) Result of simulation 

Fig. 76 shows a result of a simulation which is 
performed on the functional model of Fig. 75 and the 
respective mechanistic models using parameters of Table 13. 
Incidentally, in Table 13, an efficient 6 u of heat 
radiation of a motor coil resistance, a thermal time 
constant T M of the motor coil resistance, a coefficient 5 
P of heat radiation of a positive characteristic thermistor, 
and a thermal time constant V P of the positive 
characteristic thermistor are determined in accordance with 
an identification putting emphasis on an accuracy of a 
motor current . 
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In the top portion of Fig. 76 current I 0 _ L [A] and 
voltage V P [V] of a positive characteristic thermistor are 
shown, wherein measurement data are denoted as fine, thick 
line and calculation results are denoted as thick, weak 
line, and further a calculation result of a consumed power 
W P [W] is added. The second portion shows the measurement 
value of the motor-lock current I G _ L [A] and the current 
deviation A I ERR [A] of the calculation result. The third 
portion is the calculation result of the temperature t p 

[°C] and the resistance R P [Q] of the positive 
characteristic thermistor. The fourth portion is the 
calculation result of coil temperature t M [°C] , coil 
resistance R M [Q] , and the consumed power W M of the 
motor. Regarding the coil resistance R M [Q] , it is 
represented by a difference from the resistance R MS [Q] at 
the reference temperature 2 5 [°C] . in Fig. 76, a 
vertical dotted line denotes a resistance critical point of 
t p = 1 0 0 [°C] , and at this point or so the respective 
characteristics change suddenly. From this result, it 
would be understood that the identification results of the 
non- linear model are well coincident with the measurement 
data. 

Further to the result of Fig. 76, as seen 
therefrom, the behavior of the positive characteristic 
thermistor reaches the resistance critical point at 14.5 
[sec] , and from this point the current suddenly goes 
down, while the voltage goes up. This means that the 
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overcurrent of the motor starts to be blocked at the 
boundary of the resistance critical point. At that time, 
the block state (values of the right side) offers current 
I o _ L =0. 2 1 [A] , voltage V P =9. 6 5 [V] , and 
5 resistance R P =4 5. 9 [Q] , and this increases about 92 
times . With respect to the consumed power W P , while the 
current decreases, the voltage is high. Thus, a power W P 
= 2. 14 [W] is supplied to maintain the temperature t P — 
13 4. 2 [°C] of the positive characteristic thermistor. 
10 With respect to the state of the coil, it is understood 
that the highest temperature of the coil is t M = 7 1. 7 

[°C] and in the block area of the right side with respect 
to the resistance critical point it is lowered to the 
temperature t M =3 2. 7 [°C] and consumed power W M = 0 . .14 4 
15 [W] , so that burning of the coil is prevented. 

Identification of coefficient of heat radiation 
thermal time constant r M> coefficient of heat 
radiation 6 P and thermal time constant X P performed in a 
similar fashion to that of the above-mentioned motor 
20 characteristic. This identification is outlined as follows. 

The reason why the lock current I 0 _ L in Fig. 76 is 
gradually decreased after the compulsory stop is owing to 
an influence of coefficient 6 M of heat radiation and 
thermal time constant X u of the motor coil. Time up to 
2 5 the resistance critical point is strongly dominated by 

coefficient d P of heat radiation and thermal time constant 
X P of the positive characteristic thermistor. Times up to 
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the resistance critical point, where both are varied by ± 
1 0 [%] , are given with <5 P : 1 5 . 8 (+9. 2%) -13. 1 
(-9. 4%) [sec], r P : 1 6. 24 (+12. 3%) -12. 8 
(—11. 5%) [sec] . By way of example, in case of the 
5 motor coil, <S M is 14. 36 (-0. 7%) -14. 56 (+0. 

7%) , and T M is 14. 16 (-2. 1%) -14. 84 (2. 6%) 
[sec] . From this, it is understood that the influence 
on the resistance critical point is little. Incidentally, 
( ) is % value for time 14. 5 [sec] . 
10 The identification here performed implies that for 

the identification of the above-mentioned functional model 
the identification of the mechanistic model for determining 
the internal characteristic is performed. This method 
implies it is possible to identify individually the 
15 functional and mechanistic models including the non-linear 
elements in form of independent model. 
- - (Processing procedure of a functional model 

including the non-linearity) 

The above-mentioned contents are rearranged, and 
20 there will be described simply a basic procedure to execute 
a simulation by a functional model and a mechanistic model 
incorporated into the functional model. 

First, a non-linear parameter, incorporated into a 
linear functional model (a government equipment) and a 
25 mechanistic model (a mathematical model), are previously 
calculated and the result is substituted and renewed. 
Consequently, regarding the processing procedure, first. 
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the mechanistic model is executed to determine a 
substitution quantity, and then the linear functional model 
is executed. Fig. 7 7 shows an execution sequence of the 
functional and mechanistic models . 
5 Fig. 7 7 shows a flowchart of a basic process in 

the event that a government equation of a small scale of 
functional model is given in the discrete form to be 
executed. In Fig. 77, a repetition of the process from the 
computation start to the end is performed in unit of a 

10 sampling period in which 1 a prediction computation for an 
estimated observation quantity. As an actually used 
computer software, a whole system of functional and 
mechanistic models module is incorporated into a software 
module of the standard execution environment on a 

15 telescopic system basis. At that time, regarding the 

respective mechanistic models, modules for the mechanistic 
models prepared in form of the standard function are 
selectively called in accordance with an object of the use 
to be executed. 

20 In this flowchart, first, after execution of the 

initialization as a preparation of execution of a 
computation, the computation is started (step a). 

According to the computation, the computation is 
performed on the mechanistic model for the non-linear 

25 parameter and the switch element for the subsequent 

sampling time. Specifically, a computation is carried out 
on the mechanistic model for derivation of the non- linear 
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parameter (step b), a condition decision is carried out 
(step c), and the initialization before execution of the 
linear functional model, such as the storage element, is 
performed (step d) . 
5 Next, the non- linear parameter and the state of 

the switch element, which are determined in accordance with 
the above-mentioned manner, are substituted for the 
functional model (step e). With regard to the subsequent 
sampling time, the functional model is executed (step f ) . 

10 In the event that the subsequent sampling period 

is further continued (step g) , the process goes to a step h 
in which a prediction computation for an estimated 
observation quantity of the subsequent sampling period is 
performed in accordance with the result of computation of 

15 the sampling period now obtained, and the newly determined 
estimated observation quantity is transferred to the 
mechanistic model. 

In the mechanistic model, a computation for the 
non -linear parameter and the like related to the sampling 

20 time associated with the estimated observation quantity is 
carried out in accordance with the estimated observation 
quantity, and a result of computation thus obtained is 
transferred to the functional model. 

Repetition of those processes make it possible to 

25 perform a reproduction (a simulation) of the system 
including a non-linearity. 
(Reference) 
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1 . Symbols for modeling 

(1) Symbols for a linear model 

Regarding a block diagram representative of a 
model, a rule of the general block diagram is adopted, and 
5 a modeling is implemented by adding symbols shown in Fig. 

78. Symbols shown in Fig. 78 are used mainly for a linear 
model . 

(2) Symbols for a non- linear model 

For modeling of a non-linearity of various types, 
10 there are known the associated various expressions. In 
effect, any expression that rule and way of establishing 
modeling for content, nature and operation of a non- 
linearity are clarified to be understood on a visual basis, 
and it is possible to be transformed to a mathematical 
15 model satisfies an object of the modeling. Fig. 79 shows 
typical symbols used for modeling of a non-linearity. 

(3) Signals for a non-linearity 
(1) Quantity of operation 

Quantity of operation denotes transformed physical 
20 quantities such as absolute value, sign and square of state 
quantity other than the potential and the flow quantities; 
and logical signals such as 0 (False) and 1 (True) 
operating switch operators and the like. However, in the 
event that the switch element of the functional model is 
25 operated with the associated mechanistic model having a 
nesting structure, it is preferable to represent it with 
the subsequent substitution quantity. In some case, it 
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happens that position and displacement wherein velocity is 
integrated are also represented by the same. 
(2) Substitution quantity 

Substitution quantity denotes a line to substitute 
5 a value for a parameter from the mechanistic model. Name 
of a signal line may adopt the same name as the variable 
name for the non- linear parameter. The substitution 
quantity makes it possible that the mechanistic model is 
independent of the functional model. Selection of the 
10 substitution quantity permits the non-linear parameter and 
the mechanistic model to be telescopic system. 
(4) Logical operator 

A logical operator is used for a switch element, 
and condition decision and estimated observation quantity 

15 are used for this operation. 

(a) Estimated observation quantity 
- - There is a symbol for observing a state quantity 

of a sampling period preceding by one for a condition 
decision and the like. Regarding the state quantity 

20 designated with such a symbol, after execution of the 

government equation of the present time (k) , the estimated 
observation equation is carried out again through 
substitution of a dependent variable to an independent 
variable of the internal state quantity, so that the state 

25 quantity of the subsequent sampling time (k + 1) is 

estimated. This estimated observation quantity is one of 
concepts necessary for implementing a modeling of a non- 
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linear element in accordance with a modeling scheme now 
proposed . 

(b) Logical operator 

A logical operator used for the switch element 
will be described in detail later. Here, a switch operator 
is referred to as the switch element . 

(5) Non-linear operator 

A non-linear operator is a symbol for performing a 
modeling of a mechanistic model which is to be incorporated 
into a functional model, 1 and is used mainly for 
transformation of the observation quantity or the operating 
quantity. 

(a) Absolute value 

The absolute value is derived, by excepting sign 
of the state quantity and the operating quantity. In the 
mathematical model, it is described in expression of 
blabs (A) or the like. 

(b) Sign 

Sign of positive or negative (+) is derived from 
the state quantity and the operating quantity. In the 
mathematical model, it is described in expression of b = 
sign (a) or the like. 

(c) Square 

An operating quantity of a value of the state 
quantity and the operating quantity squared is derived. 
This symbol is the same as the symbol for the 
multiplication of two or more inputs. 



(d) Initialization (integration) 

In the initialization, a value of the integration 
symbol is initialized through a condition decision or a 
switch element. The initialization is performed before 
execution of the mechanistic model and switch element. 
However, in some case, it happens that the initialization 
is performed after execution of the mechanistic model and 
switch element. Fig. 80 shows a typical way of 
initialization. In Fig. 80, when condition of part (a) of 
Fig. 80 is not applied, the integration value x is reset to 
zero; when condition of part (b) of Fig. 80 is applied, the 
integration value x is initialized to B; and when condition 
of part (c) of Fig. 80 is not applied, the integration 
value x is initialized to B. Those mathematical models are 
expressed as follows . 

if(A*0) then ( ) else (x = 0) 

if (A = 0) then {S WA =l,x =B) 

if (A * 0) then (S WA = 1) else (S WA = 0,x = B) 

(110) 

In the model for initialization, it is inhibited 
that a model, such as an addition symbol which cannot be 
operated in the condition decision, is incorporated between 
the integration symbol and the condition decision of the 
part (a) or the switch symbols of the parts (b) and (c). 

(e) Condition decision 
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Condition decision is an observation quantity at 
the subsequent sampling time, which is observed with the 
symbol of the estimated observation quantity shown in Fig. 
79, and a symbol for generating an operational quantity of 
5 the switch element from the operational quantity generated 
within the mechanistic model. For the condition decision 
equation, there is used an operator dealing with relations 
such as large, small, equivalent etc. and executing logical 
operations. As a result of the operation, when it is 

10 applied, the operational' quantity of 1 (true) is 

outputted, and when it is not applied, the operational 
quantity of 0 (false) is outputted. In this case, 
operational variable names are applied for outline arrows 
on a colored background indicative of the operational 

15 quantity or the substitution quantity. Here, the 

mathematical model of the condition decision is described 
in the IF sentence such as C- language, which is generally 
used. 

(6) Functional equation 

20 Functional equation serves to perform a 

designation to derive functional values of the state 
quantity and the operational quantity in form of the 
operational quantity. For example, a designation that A is 
inputted into the parenthesis of the functional name s i n 

25 0 and B is outputted implies a mathematical model of B = 

s i n (A) . 

( 7 ) Operation 
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A signal to operate a model from the exterior is' 
generated. For example, there is an operation of a switch. 
(8) ON-OFF switch 

Fig. 81 shows a general ON-OFF switch in 
association with a contact type of switch which is used in 
the electric system. 

In Fig. 81, an NO switch (Normal Open Switch) 
shown in part (a) of Fig. 81 is a switch of normally 
turning off, and an N C switch (Normal close Switch) shown 
in part (b) of Fig. 81 is a switch of normally turning on. 
The mathematical model of the NO switch element shown in 
part (a) of Fig. 81 is expressed by the following equation. 

f(O p ) then {S WON = 0) else (S WON = 1)1 

(111) 

" - The mathematical model of the NO switch element 

shown in part (b) of Fig. 81 is expressed by the following 
equation . 



f(0 F ) then (S WOFF = 1) else (S W0FF = 




(112) 

The upper sides of eqs . (Ill) and (112) denote the 
condition decision equations, and the lower sides denote 
the equations of the switch element for operating the state 
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quantity. The NO switch of eq. (Ill) turns on, when the 
condition decision equation O p is applied so that the 
switch element S WON is 1 , and turns off, when the 
condition decision equation O p is not applied so that the 
5 switch element S WON is 0. On the other hand, the NC 

switch of the equation (112) turns on, when the condition 
decision equation O p is not applied so that the switch 
element S WO FF is 1 , and turns off, when the condition 
decision equation O p is applied so that the switch element 

10 Swon is 0. 

( 9 ) Logical product and logical sum 
Logical product (AND) and Logical sum (OR) are a 
logical operation associated with the logical product and 
the logical sum of Mini-Max method for the fuzzy operation. 

15 The logical product denotes a symbol having a function of 
selecting and outputting the minimum value of the entered 
state quantity. Figs. 82 and 83 show a logical product and 
a logical sum in comparison with the switch elements, 
respectively . 

20 In Fig. 82, the inside of the logical product has 

NO switch elements common in their outputs and a condition 
decision to select the minimum state quantity. In Fig. 83, 
the logical sum has a condition decision to select the 
maximum state quantity. Both the logical product and the 

25 logical sum are the same model except for the decision 
equation for the condition decision. 

The mathematical model of the logical product of 
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Fig. 82(b) is expressed by the following equation. 



*f(I A < !b) ^en (S w _ 0 = 0, S w _, = 1) else (S WJ) = 1, S w _, = 0) 

(113) 

The mathematical model of the logical sum of Fig. 
83(b) is expressed by the following equation. 

if {I a > I B ) then OVo = °> S w _ x = 1) else (S w _ 0 = 1, S w _, = 0) 

(114) 

2 . Basic function element of linear model 
(1) State quantity 

The state quantity is divided into the subsequent 
flow and potential quantities. 

- - CD The potential quantity indicates a state 

quantity indicative of a translation quantity of medium 
carrying energy such as voltage, velocity and rate of flow. 

(D The flow indicates a state quantity indicative 
of strength and quantity of energy of medium in unit 
quantity such as current, power, fluid pressure. 

The flow and the potential quantity are used in 
their combination to form energy. For example, it is 
utilized for modeling that the product of voltage and 
current is an electric power, and the product of velocity 
and power is a working factor. This implies that the base 
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of the functional model is an energy principle, and it is 
possible that the model of the system is represented by the 
potential quantity and the flow quantity, so that all the 
physical laws related to the potential quantity and the 
flow quantity are represented on a model. 

In the potential quantity and the flow quantity, 
the integrated value is expressed in form of the storage 
quantity of amount and strength of energy. The potential 
quantity and the flow quantity correspond to, for example, 
translation position and power product of the mechanical 
system. Further, as the internal energy constituting 
entropy of the thermodynamics, the temperature related to 
thermal energy for compression of gas and resistance loss 
is also the storage quantity. Here, those are referred to 
as a storage state quantity in a general term. A system 
connected with the potential quantity and the storage 
quantity is referred to as a potential system, and likely a 
system for the flow quantity is referred to as a flow 
system. Both the systems are in a duality. Those systems 
form the base of the mechanical and the electric 
engineerings. According to the mechanical engineering, 
velocity of the potential system or distance of the 
associated integration quantity (position) are dealt with 
as a result in accordance with the force (pressure) of the 
flow system. On the other hand, according to the electric 
engineering, the current of the flow system is dealt with 
as a result in accordance with the voltage of the potential 



system. According to the modeling scheme here described, 
the potential system and the flow system are dealt with on 
an equivalent basis and the duality between the potential 
and the flow systems is utilized, so that theories and laws 
of the fields of the respective engineering are mutually 
accepted to remove the barrier. 

( 2 ) Parameter 

(a) Characteristic 

Characteristic is to associate the flow quantity 
with the potential quantity and represents properties of a 
system. The characteristic is divided into two sorts of a 
storage characteristic of storing an energy and a loss 
characteristic involving an energy loss such as an electric 
resistance and a viscosity damping. Such a characteristic 
generates an energy (motion energy) represented by product 
of the potential and the flow quantities, and also 
generates an energy for heat and position generated inside. 
The internal energy is translated into the temperature and 
the distance so that it may be expressed in form of the 
output state quantity. Incidentally, the characteristic 
placed in the multidimensional space has coordinate 
information designating an operational direction inside. 

CD The storage characteristic is classified into 
two types in one of which change of the flow quantity is 
stored in form of an energy, as in the inductance and 
stiffness of a spring, and is changed into the potential 
quantity, and in another of which change of the potential 
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quantity is stored in form of an energy, as in the electric 
capacity and mass, and is changed into the flow quantity. 

(D The loss characteristic is classified into two 
types in one of which loss depending on the potential 
5 quantity, as in the leakage resistance of the electric 

capacity and the viscosity resistance by fluid, is changed 
into the flow quantity, and in another of which loss 
depending on the flow quantity, as in the electric 
resistance and the internal attenuation due to the material 
10 compression, is changed into the potential quantity, 
(b) Coefficient 

Coefficient denotes numeral and quantity giving 
association between flow quantity to flow quantity or 
potential quantity to potential quantity. The factor is 

15 classified into the following four types. 

(D Non-unit quantity associating the same physical 
unit systems with one another, such as the ratio of winding 
number of coils of an electric transformer and the ratio of 
gear wheels. In case of the factor that the value 

20 appearing on the mathematical model is 1, it means a direct 
coupling between the elements. 

(D Unit quantity associating the mutually 
different physical unit systems, such as torque constant 
and velocity constant of a motor, radius of a tire, and 

25 cross-sectional area of a piston. In the event that 

mutually different physical unit systems such as rotation 
and translational motion perpendicularly intersects one 
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another, it is necessary for the physical factor to have a 
function of coordinate transformation. 

(D Non-unit quantity for transforming a coordinate 
system of the state quantity of a link mechanism moving in 
a multidimensional coordinate space and a motion of a car. 

@ Non-unit quantity involving increment and 
decrement of an energy in which a factor is multiplied by 
one side of a pair of potential and flow quantities, such 
as a loss factor and a gain of a control system. 

With respect to the factors of the above-mentioned 
®~(3), the same factor is multiplied by both the potential 
and the flow quantities of a pair, and the transformation 
is performed in such a manner that no increment and 
decrement of energy is involved. The concept of a 
transformation factor for an application of a factor having 
meaning common to the potential and the flow quantities is 
an aspect of the functional model. This concept is one of 
reasons that make it possible to implement modeling 
exceeding the physical unit system and to permit the 
mutually different physical unit systems to exist in the 
same model. The factor serves as a weight factor effecting 
on a square basis for the storage characteristic an the 
loss characteristic. 

(c) Switch element 

A switch element is a logical element for 
connecting and disconnecting the state quantity inside the 
functional model. The functional model into which the 
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switch element is incorporated is a non-linear model for 
varying a model structure. It is considered that a factor 
value of the switch element is a non-unit quantity of 1 
(true) and 0 (false) or 1 • 0 * — 1, and is a special use of 
factor for providing the state quantity on a discrete basis, 
(d) Attached load 

An attached load serves as generation and 
absorption of the state quantity inside the mode. The side 
load is concerned to two types associated with the 
potential and the flow quantities. The side load has such 
a meaning that all single potential or flow quantities of 
no pair are regarded as the side load, and the state 
quantity not restricted is excluded form the inside of the 
system. 

CD One corresponding to the potential quantity 
such as an internal voltage of a battery, a voltage drop of 
a ^diode, and an ambient temperature of an equipment. 

(D One corresponding to the flow quantity such as 
a frictional torque, a transmission torque of a clutch, and 
a current source for supplying a constant current. 

d) A work of the side load is divided into a 
generating source for generating the state quantity, such 
as the internal voltage of the battery, and an absorption 
source for absorbing the state quantity, such as a 
frictional torque. Both the generating source and the 
absorption source are determined in work of generation or 
absorption through incorporation into a model. 



® As the especial use, there is a case where 
driver's operation and instruction are indicated in form of 
a signal. 
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